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ABSTRACT 

We present numerical simulations of electron-positron pair cascades in the magne- 
tospheres of magnetic neutron stars for a wide range of surface fields [Bp = 10"^^- 
10^^ G), rotation periods (0.1-10 s), and field geometries. This has been motivated by 
the discovery in recent years of a number of radio pulsars with inferred magnetic fields 
comparable to those of magnetars. Evolving the cascade generated by a primary elec- 
tron or positron after it has been accelerated in the inner gap of the magnetosphere, 
we follow the spatial development of the cascade until the secondary photons and pairs 
leave the magnetosphere, and we obtain the pair multiplicity and the energy spectra 
of the cascade pairs and photons under various conditions. Going beyond previous 
works, which were restricted to weaker fields {B<a, few x 10^^ G), we have incorpo- 
rated in our simulations detailed treatments of physical processes that are potentially 
important (especially in the high field regime) but were either neglected or crudely 
treated before, including photon splitting with the correct selection rules for photon 
polarization modes, one-photon pair production into low Landau levels for the e^, and 
resonant inverse Compton scattering from polar cap hot spots. We find that even for 
B ^ Bq = 4 X 10^'^ G, photon splitting has a small effect on the multiplicity of the 
cascade since a majority of the photons in the cascade cannot split. One-photon decay 
into e~^e~ pairs at low-Landau levels, however, becomes the dominant pair production 
channel when B>3x 10^^ G; this tends to suppress synchrotron radiation so that the 
cascade can develop only at a larger distance from the stellar surface. Nevertheless, 
we find that the total number of pairs and their energy spectrum produced in the 
cascade depend mainly on the polar cap voltage BpP~^ , and are weakly dependent 
on Bp (and P) alone. We discuss the implications of our results for the radio pulsar 
death line and for the hard X-ray emission from magnetized neutron stars. 

Key words: radiation mechanisms: non-thermal - stars: magnetic fields - stars: 
neutron - pulsars: general. 



1 INTRODUCTION 

The pair cascade in the magnetosphere of a pulsar has 
long been considered an essential ingredient for the pul- 
sar's nonthermal emission, from radio to gamma rays (e.g. . 



Sturrocklll97"ll:lRuderman fc Sutherlandll975l:[Meir^l2004l: 



Thompson! 20041 ) . More recently it has been suggested that 
the pair cascade is also necessary for nonth ermal cmis- 
sion from magnet ars (e. g., [Bcloborodov fc Thomps on 20o3; 
lThompson|[2008al lbl: see I Woods fc Thompsonlbood for a re- 
view of magnetars). The basic pair cascade involves several 
steps: (i) acceleration of primary particles by an electric field 
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parallel to the magnetic field; (ii) gamma ray emission by the 
accelerated particles moving along the magnetic field lines 
(either by curvature radiation or inverse Compton upscat- 
tering of surface photons); (iii) field-assisted photon decay 
into electron-positron pairs as the angle between the pho- 
ton and the magnetic field line becomes sufliciently large, 
or pair production by two-photon annihilation in weak- 
field regimes; (iv) gamma ray emission by the newly-created 
particles as they lose their transverse energy through syn- 
chrotron emission; (v) further pair production and gamma 
ray emission via steps (iii) and (iv). The dense, relativis- 
tic (Lorentz factors 7 > 100) electron-positron plasma gen- 
erated by this cascade is a required input in many mod- 

I 1 1 1 I 1 

els for the p ulsar radio emission (e.g., 'Melrose 1995, '200|; 

iBeskinI 1 19991 : rMclikidze. Gil, fc Patarava .2000 : .Lvubarskvl 
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\200i . I2OO8I : iLvutikovl |2007| ). while the high-energy pho- 
tons emitted in pair cascade models can reproduce the ob- 
served pulse profiles and phase-resolved spectra of gamma- 



taken into account (see, e.g., 


Romani & Yadig:aj-oelul 


1995 


Chenff. Ruderman. & Zhand 


2OOOI: Dvks & RudakI 


200a 


Harding et al.1 20081: Bai & SDitkovskv 20091). We note in 



passing that the dense pair plasma generated by this cas- 
cade also pl ays an impor tant role in models of pulsar wind 
nebulae (see lAror3l2007l for a review). 

The behavior of the pair cascade in the superstrong field 
regime (magnetic field strengths B>Bq = 4.414 x 10^^ G) 
and its effect on emission from pulsars and magnetars is 
somewhat puzzling. For example, of the dozen-or-so ob- 
served magnetars, only two show pulsed radio emission, 
and it is of a completely different nature than the emis- 
sion from "standard" radio pulsars (e.g., the radio pulsa- 
tions are transient and appear to be correlate d with strong 
X-ray outb ursts from the magnetars; see ICamilo et al.l 
I2OO7I I2OO8I ). In contrast, several radio pulsars with in- 



tars have been discovered (e.g., Kas 


pi & McLauehlinll2005l: 


IVranevsevic. Manchester. & Melrose 


20071). Whv the stan- 



dard mechanism for pulsed radio emission turns off for mag- 
netars but not for these pulsars is unknown. 

There have been only a few publications devoted 
to numerical simulations of the pair cascade in pul- 
sar magnetospheres. For moderate-strength magnetic fields 
(B<5 X 10^^ G), si gnificant progress has been made. 
iDaughertv fc Hardind (1982) present simulations of the cas- 
cade initiated by a single electron injected from the neu- 
tron star surface, emitting photons through curvature ra- 
diation, for (polar) surface field strengths Bp up to 5 x 
10^^ G and rotation periods P = 0.033-1 s. In a later 
paper IDaughertv fc Hardind Il996l ') they consider gamma 
ray emission from the entire open-field-line region of the 
magnetosphere, using a simplified acceleration model and 
for Vela-like pulsar parameters (B^ = 3 x 10 ^'^ G and 
P = 0.089 s). ISturner. Dermer. fc Michell (| 19951 ) present a 
similar simulation to that of Daugherty fc Harding, but 
for cascades initiated by electrons upscattering photons 
through the inverse Compton process (a gain for Vela-like 
parameters). iHibschman fc Aronj (|2001bl ) develop a semi- 
analytic model of the inner gap cascade, both for curva- 
ture radiation-initiated and inverse Compton scattering- 
initi ated cascades, applicab le for -6^3 x 10^^ G (see 
also IZhang fc Hardind I2OO0I ). Cascades occurrin g in the 
outer magnetosphere have also been simulated, by Romanil 



l|l996l) for Vela- and Crab-like (B^ = 4 x 10'^ G and 
P = .0 33 s) parameters (see also Cheng. Ho. fc Rudermaiil 
ll986al lbl: ICheng. Ruderman. fc Zhandl2000l ). 

However, for superstrong magnetic fields {B > Bq = 
4.414 X 10^^ G) only limited aspects of the full casc ades have 
been studied. For example, lArendt fc Eilekl l|2002l ) simulate 
the cascade for Bp < 10^^ G and P = 0.033 s (for both a 
pure dipole and a more complex field geometry), but with 
the simplification that all photons radiat ed by the primary 
partic le are e mitted from the surface. iBaring fc Hardind 
liool) (see also lHarding. Baring, fc GonthieJll997l ) use this 
same simplification to study the effects of photon splitting on 
the cascade for field strengths up to B = 2x 10^* G (however, 
they assumed that both photon modes can split, and thus 



overestimated the e ffect o f photon splitting; see Section [3]2]). 
tearing fc Hardind (|2007l ) model the process of resonant in- 
verse Compton scattering of photons from the neutron star 
surface (with the blackbody temperature T — 6 x 10^ K) in 
the same field range, b ut only for single scattering events 
(see also iDermeij Il990l ). The magnetosphere acceleration 
zone in the superstrong, twiste d field regime of magnetars 
is inv estigated analytically by iBeloborodov fc ThompsonI 
liool) for cascades occurring i n the closed field l ine region 
of the magnetosphere and by iThompsonI (|2008al lbl) in the 
open field line region. 

In this paper we present numerical simulations of the 
pair cascade from onset to completion. Motivated by the 
lack of full cascade results for the superstrong field regime, 
and in light of the unexplained differences between the ob- 
served emission properties of high-field radio pulsars and 
magnetars, we run our simulations in magnetospheres with 
field strengths up to 10^^ G. We consider several impor- 
tant factors that affect high-field cascades, including pho- 
ton splitting, pair creation in low Landau levels, photon 
polarization modes (_L or || to the magnetic field direc- 
tion), and resonant inverse Compton scattering. We use 
our simulations to generate spectra of the high-energy pho- 
tons and the electron-positron plasma produced by the cas- 
cade. Additionally, we use our simulations to comment on 
the conditions for when the radio emission mechanism no 
longer operates in the neutron st ar magnetosphere, the so- 
called "pulsar death line" (e . g.. iRuderman fc Sutherland! 
1975!; 'Chen & Ruderman" 1993^; 'Hibschman & Arondl2001 
Harding fc Muslimov 2002; Harding, Muslimov, fc Zhan!_ 
2002! : [Medin fc LaiH2007i ). While the results of our simulation 



are most applicable to cascades occurring in the open field 
line region of the magnetosphere (since the primary particles 
are injected into the magnetosphere along open field lines), 
some of our results are also relevant to cascades occurring in 
the closed field line region for magnetars, e.g., the products 
of a cascade initiated by a photon injected into a non-dipole 
magnetosphere. 

A necessary component of any pair cascade simulation 
is a model of the magnetosphere acceleration zone, or "gap" , 
where the cascade originates. In real magnetospheres of pul- 
sars and magnetars, the acceleration of primary particles is 
coupled to the rest of the cascade (e.g., charged particles pro- 
duced in the cascade can screen out the acceleration poten- 
tial). However, there is significant uncertainty about the pre- 
cise nature of the acceleration gap. A number of models have 
been proposed for the location of the gap, from inner mag- 
netosphere accelerators (bot h "vacuum" and "space-charge- 
limited flow" types; see , e.g.. [Ruderman fc Sutherland 19751: 



Arons fc Scharlemann 19791: Muslimov fc Tsvgan 



1992 



Hibschman fc Aronsll2001al : iMedin fc Laill2007l : rrhompsonl 
2008a|[b|). t o oute r magnetosphere acceler ators (e.g. , 
Cheng et all Il986a| [bl: [Romani 1996; Cheng et al.1 I2OO0I : 



Takata et aL[|2006l ). to hybrid inner-outer magnetosph ere ac- 



celer at ors ( "slot" gaps and extended oute r gaps; e.g. , lArong 
1 19831 : [Muslimov fc Hardind |2003| . [20041 : iHirotanil I2OO6I ). 
Non-steady (oscillat ory) inner ga ps have also been dis- 
cusse d recently (e.g . . ISakai fc Shibata,200A [Levinson et al.l 
I2OO5I : iBeloborodovl [20081 : ILuo fc Melrosd [2OO8I ). Numer- 
ical simulations of force-free global magnetospheres in- 
cluding magnetic-field t wisting near the light cylinder 
have been performed (e.g.. [Contopoulos. Kazanas. fc Fendt[ 
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19991: iGruzinwl l2005l: ISpitkovskvl [ioo^: iTimokiiiiil l2006l : 



Komissarovl 20061 : Kalapotharakos fc ContopoulosI |2009^ ■ 
but they do not yet include any particle acceleration, or 
pair creation self-consistently. Therefore, in this paper we 
decouple particle acceleration from the rest of the cascade 
and focus on the cascade produced by a primary electrorQ 
injected into the magnetosphere with a Ri ven initial Lorentz 
factor 7o (cf . iDaughertv &: HardintJI 19821 ) . We also consider 
the cascade produced by a single "primary" photon emitted 
by the primary electron, in the case where photon emission 
within the acceleration gap is important (i.e., for cascades 
where the dominant mechanism for high-energy photon pro- 
duction is inverse Compton scattering). 

The outline of the paper is as follows. In Section [2] we 
summarize our method for estimating the initial parameters 
(e.g., 7o) of the primary cascade particles, for use in our sim- 
ulations. In Section [3] we describe the details of the numer- 
ical simulations, both for cascades with photon production 
dominated by curvature radiation and by resonant inverse 
Compton scattering (resonant ICS, or RICS). In Section [4] 
we present our results (e.g., photon and pair plasma spec- 
tra) for a wide range of parameters: surface magnetic fields 
B = lO^'^-lO^^ G, rotation periods 0.1-10 s, surface tem- 
peratures T = (0.3-3) X 10^ K, and pure dipole and more 
complex field geometries. In Section [5] we summarize our 
findings and discuss their implications for the radio emis- 
sion and high-energy (hard X-ray and gamma-ray) emission 
from pulsars and magnetars. Some technical details (on our 
treatment of inverse Compton scattering, on our treatment 
of attenuation coefficients and e^e~ energy levels for pair 
production, and on deriving semi-analytic fits to our nu- 
merical results) are given in the appendix. 



2 ESTIMATING THE INITIAL PARAMETERS 
FOR THE PRIMARY PARTICLES 

2.1 Primary electrons 

In our cascade simulation (described in Section |3} we do not 
include an actual acceleration region, since we wish in this 
work to remain as model-independent as possible. Instead, 
we model the effect of this region on the cascade by giving 
the primary electron an initial energy 7omeC^ equivalent to 
the energy it would reach upon traversing the entire gap, 
and injecti ng it into the magnetosphere at the neutron star 
surface (cf. IDaughertv fc Hardinelll982h . Obviously, this ap- 
proximation excludes a proper treatment of the slot gap and 
outer gap acceleration models. However, in most parts of the 
polar cap region (i.e., except for the boundary region adja- 
cent to the open field lines), the main voltage drop occurs 
near the stellar surface, regardless of the nature of the "gap" 
(vacuum gap or space charge limited flow). Such inner gap 
models and other global models with near-surface accelera- 
tion are allowed in our analysis. 

For a dipole magnetic fleld geometry, most active pul- 
sars with inner gap accelerators have gap voltage drops in 



the range "I> ~ ( 1-2) x 10^'' V, regardless of the accelera- 
tion model (e.g.. iHibschman fc Aron3l2001al : iMedin fc Lail 
I2OO7I . hereafter ML07). For the surface field strengths we 
are considering, B > 10^^ G, the primary electrons are not 
radiation-reaction limited with in these gaps (ML07; cf. the 
millis econd pulsar models of iHarding. Usov. fc Muslimovl 
I2OO5I ). so we can set 70 — e^/m^c?. We therefore restrict 
70 to the range (2-4) x 10*^ for dipole fields. Note that these 
large voltage drops do not occur in pulsars where the gap 
electric field is fully screened due to inverse Compton scat- 
tering by the primary electron. We discuss this case in Sec- 
tion [121 

The voltage drop across the gap can be no larger than 
the volta ge drop across the entire polar cap of the neutron 
star (e.g.. iRuderman fc Sutherlandlll975h : 



~2c c~ 



= 7 X 10''Bp,i2Po" V, 



(1) 



where R is the radius of the star (assumed in this pa- 
per to be 10 km), Pq is the spin period in units of 1 s, 
and _Bp = 10^^i3p,i2 G is the polar surface magnetic field 
strength. If the voltage drop, <1>, required to initiate pair 
cascades is not available, i.e., "I> > 3>cap, the magnetosphere 
should not produce pulsed radio emission; the locus of points 
where $ = "I>cap defines the pulsar death line0 A typical 
death line for an inner gap model, plotted in P-P space, is 
shown on the left panel of Fig. [T] The line was made using 
three assumptions: (i) The magnetosphere field geometry is 
dipolar, (ii) The pair cascade occurs primarily above the 
gap, (through curvature radiation) once the primary elec- 
tron has reached a large Lorentz factor 70 ~ 10^. (iii) The 
spindown power of the pulsar, given by 



E - 



-mo. - 



— p3— 



(2) 



is approximately equal to the spindown power of a magnetic 
dipole with its magnetic field and rotational axes orthogonal 
to each other: 



E'. 



3 • 

The polar magnetic field strength inferred from this 
frequently- used approximation is 



B. 



p. 12 



2.0\/PoP-i5, 



(4) 



where P-15 is the period derivative in units of 10~^^ s/s and 
/ — 10*^ g-cm^ is assumed. 

A well-known problem with the death line made using 
these assumptions is that it cuts right through the middle of 
the main grou p of p ulsars (Ruderman & Sutherland"l975|; 
Hibschman fc Arond [loOla; Harding fc Mushmov 2003; 
Medin fc Laill2007l ): i.e., the model incorrectly predicts that 
there will be no radio emission from many neutron stars that 
are observed to be active pulsars. 

Several authors have proposed models of the neu- 
tron star magnetosphere that shift the theoretical death 
line closer to the observed death line by altering one 



^ Although the primary particle could also be a positron (or even 
an ion), we assume here for simplicity that the pulsar is oriented 
such that electrons are accelerated away from the star. 



^ This applies to the vacuum gap model. In the space-charge- 
limited-flow model the condition is 3> = Kg<&cap, with Kg ~ 0.15 



(e.g.. [Hrbschman fc Aron£ir2001ar) . 
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Figure 1. Pulsar death lines. Death lines are shown for pulsars with dipole magnetic fields, dipole fields offset from the center of the star 
by Ar = 0.95-R, and magnetic fields with extended polar caps 100 times larger than the dipole value 9cap = \/flR/c (left panel); and for 
pulsars with magnetic field curvatures TJc = -R at the surface (right panel) . Note that these death lines do not apply for the millisecond pul- 
sar population in the lower left corner of the diagram, as their short periods and low magnetic field strengths cause the primary electron to 
be radiation reaction limited. In each panel, rotation-powered pulsars (ATNF catalog, http://www.atnf.csiro.au/research/pulsar/psrcat) 
are labeled by crosses, while magnetars (McGill catalog, http://www.physics.mcgill.ca/~pulsar/magnetar/main.html) arc labeled by 
solid circles and the two radio magnetars are labeled by solid triangles. 



or more of the assumptions made above. In some mod- 
els the magnetosphoro geometry is not a centered dipole, 
but instead is an ofltsot dipole (Arons 199^) or a twisted 
dipole (e.g.. [Thompson. Lvutikov. fc Kulkarnil2002l ). or con- 
tains quadrupole or higher multipole components (e.g., 
IPavan et aLll2009l ). Numerical simulations suggest that due 
to twisting of the field lines near the hght cylinder the po- 
lar cap is slightly l arger than in the pure dipole case (e.g., 
IContopoulos et al] Il999l : ISpitkovskvl |2006| ) : the polar cap 
could be significantly larger if, as the pulsar spins down, 
field line reconnection (conversion of open field lines into 
closed field lines) i s too slow to keep p ace with the expand- 
ing light cylinder (jContopoulosI I2OO5I ) . All of these models 
increase the size of the theoretically allowed P-P space for 
pulsars by decreasing the radius of curvature of the mag- 
netic field lines Tic in the region where the pair cascade 
occurs, i.e., above the polar cap. Since pair creation is more 
efficient along tightly-curved field lines, a smaller radius of 
curvature allows the cascade to occur at a lower 70. How- 
ever, the magnetosphere must be highly non-dipolar near 
the polar cap in order for the models to include all pul- 
sars on the active side of the death line. For example, the 
offset dipole model can fully match observation only if the 
dipole is offset by 0.95-R or more, while an expanded po- 
lar cap (due to a twisted field or delayed field reconnec- 
tion) must be hundreds of times larger than the pure dipole 
cap; see Fig. [T] It is uncl ear whether such a strongly non- 
dipolar field is stable (cf. iBeloborodov fc Thompson! |2007| . 
where a field twist > 1 radian is unstable even in mag- 
netars). In addition, observations of pulsar radio emission 
suggest that this emission is com ing from a purely dipol e 
region of the magnetosphere (e.g.. iRankin fc Wrightll200a ). 
so any non-dipole structure at the surface must give way to 
a dipole configuration at the altitude where radio pulses are 
generated (typica lly about 1% o f the light cylinder radius; 
e.g., .Gangadhara fc Guptal l200il ; iDvks. Rudak. fc Harding 
In our simulation, we model the effects of a non-dipole 
magnetosphere by giving the magnetic field a large curvature 



[He = R; cf. iRuderman fc Sutherlandlll975l ) at the surface 
and a dipole curvature at higher altitudes (r > 2R) ; the im- 
plementation of this non-dipole model into our simulation 
is discussed in more detail in Section [3] The death line in 
this approximation is shown on the right panel of Fig. [1] 
Because of the increased efficiency of the cascades in these 
highly-curved magnetic fields, the Lorentz factor of the pri- 
mary electron upon emerging from the gap is a factor of ~ 10 
lower than in the dipole case, in the range 70 ~ (2-4) x 10® 
(iHibschman fc Aronsll2001al ; iMedin fc Laill2007l ). 

In some models the pair cascade occurs primarily within 
the gap, due to efficient inverse Compton scattering by the 
primary electron, rather than above the gap. These cas- 
cades occur at much lower energies of the primary electron 
(7 ~ 10^-10**), such that all observed pulsars can provide 
the voltage necessary to initiate this type of cascade. How- 
ever, ICS cascades are generally very weak, producing < 10 
electron-positron pairs per primary electron, while most ra- 
dio emission models assume a secondary particle density 
100-1000 times that of the primary electron be am (e.g., fo r 
the development of a two-stream instability; see lUsovll2002l ). 
The inclusion of this type of cascade into our simulations is 
discussed below, in Section [2.21 

In some models the spindown power of the neutron star 
differs from t hat of an or thogonal magnetic dipole, Eq. @. 
For example, IContopoulo s & Spitkovsky (2006) show that 
for nearly-aligned pulsars (with the angle between the mag- 
netic field and rotation axes a<30°), the spindown power 
can be approximated by 

This equation is nearly the same as Eq. Q for young pul- 
sars, where $ <^ $cap. However, for pulsars near death, with 
$ < $cap, the spindown power is much lower for a given po- 
lar cap voltage. Conversely, for a given P and P [which de- 
termines the observed spindown power, Eq. ((2])] the polar 
cap voltage is much larger than would be assumed by us- 
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ing Eq. ((3]). Depending on alignment and how close the gap 
potential drop "I> is to ^cap, many pulsars which were pre- 
dicted to be dead may actually have a potential drop large 
enough to generate pair cascades ($ ~ 10^ V). According to 
Contopoulos & Spitkovsky, the standard death line shown 
in Fig. [T] is consistent with the observed P-P values for all 
pulsars if the magnetic inclination angles a of nearly-dead 
pulsars are weighted towards a = 0. If this is the case, we 
can use 70 — (2-4) x 10^ for all pulsars and do not need to 
invoke a strongly-curved magnetic field geometry {TZc — R) 
or efficient inverse Compton scattering by the primary elec- 
tron in order to reproduce the observed pulsar death line. 



2.2 Primary photons 

As the primary electron traverses the acceleration re- 
gion it can pass through two distinct regimes (see, e.g., 
iHibschman fc ArondbOOlal ). First, for Lorentz factors 7 ~ 
lO'^-lO*, the electron efficiently upscatters photons through 
the inverse Compton process. Second, for Lorentz factors 
7>10^, the electron efficiently emits curvature radiation. 
The above treatment of the acceleration zone (Section 12. ip 
is best suited for pulsars where the electrons reach the sec- 
ond regime. In that case we can safely ignore the contribu- 
tions to the cascade made by photons emitted before the 
primary electron reaches full energy (70), since the number 
and energy of photons emitted through curvature radiation 
increase strongly with 7 (i.e., Nck oc 7, Ecu oc 7^^). The ap- 
proximation is poor, however, if inverse Compton scattering 
and subsequent pair production within the gap is efficient 
enough to screen the accelerating potential before the elec- 
trons can reach the second regime. In that case the photons 
produced in the gap are critical to the cascade, while the 
photons produced above the gap have a negligible effect on 
the cascade (the upscattered photons must travel a finite 
distance before pair production in order to screen the gap; 
in that distance the primary electron is accelerated to above 
resonance and exits the gap with a Lorentz factor between 
the first and second regimes of efficient photon production). 

Because ICS is strongly peaked at resonance, primary 
electrons traveling through this second type of gap will emit 
a large number of photons at a characteristic "resonance" 
energy and very few at other energies. The effect of this type 
of gap on the cascade is better model ed by A^o photons of en - 
ergy eo emitted from the surface (cf. lArendt fc Eilekll2002h . 
rather than one electron of energy 7omeC^. We therefore run 
a second version of the simulation, this time tracking the 
cascade initiated by a "primary" photon. The quantitative 
results of this simulation can be multiplied by No to obtain 
the full cascade results (e.g., the number of electron-positron 
pairs produced per primary electron). 

For a primary electron resonantly upscattering primary 
photons, we estimate the value of eo as follows. When the 
primary electron reaches a Lor entz factor 7, it upscatters 
photo ns to a mean energy (e.g.. iBeloborodov fc ThompsonI 
\2OOt ) 



e = 7 1 



v/1 + 2/3q 



mn^c 



(6) 



where /Bq = B / Bq is the ratio of the magnetic field strength 
to the critical quantum field strength, Bq — 4.414 x 10^^ G. 



The primary electron is most efficient at scattering photons 
when 



7 = 7crit 



c/kT, 



(7) 



where T is the surface temperature of the star and Ec = 
heB/rrieC is the electron cyclotron energylf] Therefore, for 



typical energy of a scattered photon is 

ernes ^ TOBp.iaTg- V Wq) MeV , 



(8) 



where Te is the surface temperature in units of 10^ K and 
/ (Pq) ~ 1 — + 2/3q is evaluated at the surface. Setting 
eo = eRics and assuming a Te range of 0.3-3, we obtain eo 
in the range 1-10 MeV at _Bp,i2 = 1 up to (0.4-4) x 10^ MeV 
at Bp,i2 = 1000. 

The number of resonant ICS photons scattered by the 
primary electron is more difficult to estimate, since it de- 
pends on the acceleration mo del. Inner gaps with space- 
charge-limited flows have (e.g.. IHibschman fc Aroii3l2001al : 
[Mcdin 2008. ') 



No 



WB-i^Po 



3/4y5/2 



(9) 



Inner vacuum gaps, with accelerating electric flelds on the 
order of 6 times larger, have No values at least 20-100 times 
smaller (the primary electron is more rapidly accelerated out 
of resonance; see, e.g., ML07). 

Note that we can also use this second cascade simula- 
tion as a diagnostic tool for the main simulation. For exam- 
ple, we can study the partial cascade initiated by a single 
curvature radiation photon emitted at some altitude in the 
magnetosphere to understand how the strength of the local 
magnetic field affects the cascade. The characteristic energy 
of curvature photons is 



(10) 



For dipole fields the typical curvature photon has an energy 
ecR< 10^-10" MeV (for 7 < 70 ~ lO'^'^), while for fields 
with TZc = R we have ecR < 10^-10^ MeV (for 7 < 70 ~ 
lO^-'). 



3 NUMERICAL SIMULATION OF PAIR 

CASCADES: PHYSICS INGREDIENTS AND 
METHODS 

The general picture of the pair cascade as modeled by our 
numerical simulation is sketched in Fig. (2] At the start 
of the simulation, an electron with initial Lorentz factor 
70 ~ lO'^'^-lO'^-^ (Section[2TTJ travels outward from the stel- 
lar surface along the last open field line. As it travels it emits 
high-energy photons through curvature radiation or inverse 
Compton upscattering. The simulation tracks these photons 

^ Note that the actual resonance condition is ei7(l — fUcosifi) = 
ec, where e^ ~ kT is the initial (before scattering) photon energy, 
/3 = \/l — 1/7^ is the ratio of the electron speed to the speed 
of light and ip is the incident angle of the photon with respect 
to the electron's trajectory. However, because the scattering rate 
depends inversely on 7 (see Appendix|A]l, photons with cos 1/) <^ 1 
are far more likely to scatter off the electron than photons with 
cos i/) ^ 1 ■ 
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NS surface ' ^^"^ 
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Figure 2. A schematic diagram sliowing the magnetospherc pair 
cascade, from initiation by a high-energy electron to completion. 
Photon splitting is also shown. The inset shows the beginning of 
a cascade initiated by a photon upscattered through the inverse 
Compton process. 



as they propagate from the point of emission through the 
magnetosphere, until they decay into electron-positron pairs 
through magnetic pair production or escape to infinity. In 
the superstrong field regime, the photon (if it has the cor- 
rect polarization; see Section 13. 2p also has a finite proba- 
bility of splitting into two photons before pair production, 
in which case we follow the two child photons in a similar 
way. The electrons and positrons created by these photons 
are tracked as they radiate away their transverse momenta 
through synchrotron radiation and then gradually lose their 
forward momenta through inverse Compton scattering. Sub- 
sequent generations of photons and electrons /positrons are 
also tracked, in a recursive manner, and the total numbers 
and energies of photons and electrons -I- positrons that es- 
cape the magnetosphere are recorded. We track each parti- 
cle until it is destroyed or reaches a height comparable to 
the hght cylinder radius. Tracking particles out to the Ught 
cylinder is most important for the primary electron, since 
although in general there is no discernible pair production 
above the radius r ~ lOR, curvature radiation continues 
up to very high altitudes (albeit very weakly, with photon 
energies 10 MeV near the light cylinder). 

In the second version of our cascade simulation, a pho- 
ton with energy eo ~ 10^-10^ MeV (see Section [2^ travels 
outward from the last open field line at the stellar surface, 
in the direction tangent to the magnetic field at that point. 
While in theory the photon should be emitted at an an- 
gle to the field line of A0ph<l/7e, where 7e is the final 
Lorentz factor of the electron after emitting the photon, in 
practice this does not matter, as there is no change in the 
final products of the cascade whether we use ASph = l/"/e 
or AOph = 0. (This is true even for photons upscattered 
by resonant ICS at B > Bq; see Section |3]3]2]). The primary 
photon and subsequent generations of electrons/positrons 
and photons are tracked in the same way as in the main 
cascade simulation. The cascade as modeled by the second 
simulation is sketched in the inset of Fig. (2] 



The input parameters for our simulation are the initial 
energy of the electron (7omeC^) or photon (eo), its initial po- 
sition (in most cases, the intersection of the last open field 
line with the stellar surface), the general pulsar parameters 
(surface magnetic field strength Bp — 10^^-10^^ G, rotation 
period P — 0.33-5 s, and surface temperature T — 10^ K 
or 5 X lO'' K), and the geometry of the magnetic field. In 
each run of the simulation, the magnetic field structure is 
given by one of two topologies: (i) a pure dipole field ge- 
ometry; or (ii) a more complex field geometry near the stel- 
lar surface which gradually reverts to dipole at higher alti- 
tudes (a non-dipole, or "multipole" field geometry). Mod- 
eling the dipole field geometry is s traightforward (see, e.g., 
Ivon Hoensbroech. Lesch. fc Kunzm998 ). but there is no ob- 
viously correct way to model the geometry for the multipole 
field case (see Section [2. Two features of a multipole field 
geometry have a strong efi^ect on the pair cascade dynam- 
ics and must be incorporated into our model: First, the ra- 
dius of curvature TZc is much smaller than dipole (we choose 
TZc = R, the stellar radius) near the surface of the star. This 
leads to a much larger number and peak energy of photons 
emitted through curvature radiation than in the dipole field 
case. Second, as a photon propagates through the magne- 
tosphere the angle between the photon and the field, which 
scales like AQph ~ Sph/TZc, where Sph is the distance trav- 
eled by the photon from the point of emission, grows much 
faster than dipole. This leads to a much more rapid decay of 
photons into pairs than in the dipole case. The integration of 
these two features into our model is discussed in the relevant 
subsectio ns below (Sectionl3.1lan d Section|321 respectively). 
Note that lArendt fc Eilekl (|2002l ') consider the first aspect of 
a multipole field geometry in their model (that TZc = R) but 
ignore the second. In all of the simulation runs we assume 
that the local magnetic field strength varies as in the dipole 
case, 



Bir,e,(t>)^Bp - 

1 7- 



rV Vsl 



' + 1 



(11) 



where (r, 9, (p) are the spherical coordinates (with the mag- 
netic north pole at r = _R and 9 = 0). Our approximation 
therefore ignores any amplification of the field strength near 
the surface caused by the complex topology. 

For simplicity we consider a "two-dimensional" cascade 
model in which all photons are emitted and travel in the 
plane defined by the local magnetic field line. Both the pho- 
tons and the electrons/positrons are tracked in the "coro- 
tating" frame (the frame rotating with the star), and any 
bending of the photon path due to rotation is ignored - this 
is expected to be valid since the cascade takes place far in- 
side the hght cylinder. Thus we shall also call this corotating 
frame the "lab" frame for the remainder of the paper. With 
this approximation the particle positions and trajectories 
are defined only in terms of r and 6 in our simulation. We 
justify this approximation below (Sections 13. II and 13. 2[1 . As 
an additional simplification we ignore any effects of general 
relativity on the photon/particle trajectory. 

The cascade simulation can naturally be divided into 
three parts: (i) the propagation and photon emission of the 
primary electron; (ii) photon propagation, pair production, 
and splitting; and (iii) the propagation and photon emission 
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of the secondarjifl electrons and positrons. Each of these as- 
pects of the simulation is described in a separate subsection 
below. At the end of this section, cascades initiated by pri- 
mary photons are discussed. 



3.1 Propagation and photon emission of the 
primary electron 

In our cascade simulation, the primary electron starts at the 
position {ro,Oo) = {R,9o) (i-e., at some angle 9o from the 
magnetic pole on the neutron star surface) with the initial 
energy ■yomeC^, and moves outward along the local mag- 
netic field line. The initial position of the primary electron 
is chosen so that it moves along the last open field line, 
whose location at the surface is given by the polar cap an- 
gle: Oo = Scap = ^/R/r-Lc, where tlc = c/Q, is the light 
cylinder radiuslf] 

The primary electron moves outward along the field line 
in a stepwise fashion. The lengths of the steps As(r) are 
chosen so that a uniform amount of energy A7 (we choose 
~ O.OOI70) is lost by the electron in each step (7 — >■ 7 — A7): 



As(r-) 



A7 



/ds 

For an electron emitting curvature radiation, 



d7 



ds 3'^ 



2 4Q/ao 



(12) 



(13) 



where o?/ = /{he) is the fine structure constant and ao is 
the Bohr radius. For a dipole field the radius of curvature is 
given by 



7^c = 



1 + 3cos^ 



\3/2 



sin 9 3-1-3 cos^ 9 



(14) 



while for a near-surface multipole field we use TZc ~ R- As 
discussed in Section [2. II we do not consider photon emission 
due to inverse Compton scattering here, since this process 
is very inefficient once the primary electron has reached the 
energy 7omeC^. We do, however, consider in our simulation 
the photon emission due to ICS by the secondary electrons 
and positrons (see Section I3.3|l which typically have 7 ^ 
70. Note that we also indirectly include ICS in our second 
cascade simulation (described in Section [3. 4p . which models 
photon-initiated cascades, by choosing photon energies eo 
that are typical of ICS photons. 

As the electron moves a distance As along the field it 
emits photons with energies divided into discrete bins (our 
simulation uses ~ 50 bins). The energy in each bin, e, is a 
constant multiple of the characteristic energy of curvature 

■* In this paper we use the term "secondary" to refer to any cas- 
cade particle except the primary electron, positron, or photon 
that initiates the cascade. The fourth generation of electrons and 
positrons, e.g., are all "secondary" particles. 

^ Since the last open field line is also the most-tightly curved, 
this choice for the primary electron's location gives us the largest 
possible cascade in our simplified model. However, in a real mag- 
netosphere there is no cascade at all along the last open field line 
and weak cascades for field lines very near the last (6 < ©cap ) , 
since the potential drop is zero on the boundary of the open field 
region. The strongest cascades occurs on field lines neither at the 
edge of the open field region nor at the center (where TZc — )■ 00). 



photons ecR ~ S^^hc/{2TZc), with e/ecu in the range 10""*- 
10. The number of photons in a given energy bin emitted 
in one step is given by the c lassical spectrum of curvature 
radiation ("e.g.. |jacksonlll99^ . 



AN, 



^ dN ^ y/SofAs -yAe^ f ^ \ 
de ~ 2tt TZc e V '^CR / 



(15) 



where Ae is the spacing between energy bins and the values 
of TZc and 7 used are averages over the interval As. Here, 
F{x) — X K^/-j,(t)dt and ^"5/3(0::) is the n = 5/3 Bessel 
function of the second kind. Note that F{x) oc x^^ ^ for x <^ 
1, and F(x) oc ^e'"" for x > 1 (e.g.. lErbeij|l966l ). 

The photons are emitted in the direction nearly tan- 
gent to the field line at the current location of the electron 
(r, 9). For a dipole field geometry the angle between the local 
magnetic field and the magnetic dipole axis is given by 



x{9) = 9 arctan 



/ tan 6* 



(16) 



see Fig. O There is an additional contribution to the emis- 
sion angle of ~ I/7, due to relativistic beaming. In reality 
this beaming angle is in a random direction; however, for our 
two-dimensional approximation it can only be in the plane 
of the magnetic field. The photon emission angle is given by 
the (projected) sum of these two angles: 



Oph = X + - cos n , 

7 



(17) 



where 11 is a random angle between and 27r. Note that ig- 
noring the three-dimensional aspect of the photon emission 
introduces an error in the emission angle of order I/7. This 
affects the location at which the photon decays (into pairs) 
in our simulation, since photon decay depends strongly on 
the intersection angle between the photon and the mag- 
netic field (see Section [3.21 below). However, as the photon 
propagates through the magnetosphere these errors (which 
are on the order of I/7 ~ 10~^ for curvature photons and 
10"'' for resonant ICS photons) quickly become negligible in 
comparison to the photon-magnetic field intersection angle, 
which grows like Sph/TZc (and so reaches the angle I/7 by 
Sph ~ 10~^R for curvature radiation and ~ O.IR for RICS). 

We also use Eq. (|17|l for simulation runs with a multi- 
pole field geometry. This is obviously a simplification, but 
we have found that in practice the photon propagation direc- 
tion has little effect on the overall cascade product (as long 
as it points generally outward). Far more important for the 
cascade is how the angle between the photon and the mag- 
netic field changes as the photon travels. As is discussed in 
Section 13.21 we artificially force this angle to change more 
rapidly with distance than in the dipole case, to account for 
the effect of the stronger field line curvature. 

The total energy lost over each step is 



eAN, ~ A7meC^ 



(18) 



Only one photon is tracked for each energy bin e at each step 
As, so the photon is given a weighting factor AA^^. In addi- 
tion to its initial position (the position of the electron at the 
point of emission r,9) and propagation direction (Oph), the 
photon has a polarization direction. For curvature radiation 
the polarization fraction is between 50% and 100% polarized 
parallel to the magnetic field curvatur e, depending on pho- 
ton frequency (|jacksonl 119981 ; see also iRvbicki fc LightmanI 
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Figure 3. A schematic diagram illustrating the photon emission 
angle. The direction of the magnetic dipole axis is given by fi. 
The electron (positron) follows the curved field line to the point 
(r, 9), then emits a photon in a cone of width I/7, inclined with 
respect to the magnetic axis by an angle x- 

1 19791 ). We therefore randomly assign the photon a polariza- 
tion in the ratio of one _L (perpendicular to the field) to 
every seven || (parallel to the field) photons, corresponding 
to 75% averaged parallel polarization. 

3.2 Photon propagation, pair production, and 
splitting 

In our simulation, the photon is emitted/scattered from the 
point (ro,ph, 6'o,ph) with energy e, polarization || or _L, and 
weighting factor AA^e (to represent multiple photons; Sec- 
tion [STT])- It has an optical depth to pair production, r, and 
to photon splitting, Tgp, both of which are set to zero at the 
moment of the photon's creation. The photon propagates in 
a straight line from the point of emission, at an angle Oph 
with respect to the magnetic dipole axis. Note that in the 
corotating frame (which is the frame we are working in for 
most of our simulation; but see Section [3. 3|) the path of the 
photon is in reality curved, with the angular deviation from 
a st raight line growing approximately as S ph^/c = Sph/rLC 
(cf. iHarding. Tademaru. &: Espositdfl978l '). Like the beam- 
ing angle (Section l3.ip . this curved path modifies the growth 
of the photon-magnetic field intersection angle and the loca- 
tion of photon decay in our simulation. However, the total 
intersection angle grows much faster with photon distance 
Sph than the deviation does (~ Sp^/TZc versus Sph/rLC, or a 
factor of r-Lc/Tlc — lOOPg''^ larger for dipole fields), so we 
can safely ignore this deviation. 

In each step the photon travels a short distance through 
the magnetosphere, Asph < 0.05rph, where (rph,0ph) refers 
to the current position of the photon; our method for choos- 
ing the value of Asph for a given photon is discussed at 



Figure 4. A schematic diagram for deriving the angle between 
the photon and the magnetic field, ip. The direction of the mag- 
netic dipole axis is given by ft. The photon propagates through 
the magnetosphere with angle 0ph with respect to the magnetic 
axis [see Eq. I I17II ] . The local magnetic field makes an angle x with 
respect to the magnetic axis [Eq. I|16|l ]. 

the end of this section. At the new position the change in 
the optical depth for pair production, Ar, and for photon 
splitting, Atsp, are calculated: 

Ar~ Asph7?||,x, (19) 

Ar,p ~ Asphi?||Pj^ , (20) 

where -R||,x = -R(| xsini/) is the attenuation coefficient for 
the II or _L polarized photons, ^}} is the angle of intersection 
between the photon and the local magnetic field, and R' 
is the attenuation coefficient in the "perpendicular" frame 
(i.e., the frame where the photon propagates perpendicular 
to the local magnetic field). For a dipole field geometry the 
intersection angle is given by 

^ = x(V) - eph , (21) 

where Oph is given by Eq. (|17p and x(^ph) is the angle be- 
tween the magnetic axis and the magnetic field at the cur- 
rent location of the photon [Eq. (|16|) ]: see Fig.|4]for a sketch. 
For the near-surface multipole field geometry we set 

tan7/-=^ = ^. (22) 
Tic R 

This approximation has the advantage of accounting for the 
effect of a strong field curvature on the photon propagation 
without requiring knowledge of the actual field topology. 

The total attenuation coefficient (in the perpendicular 
frame) for pair production is given by (suppressing the sub- 
scripts II, -L) R' — TijkRjk, where R'j^. is the attenuation 
coefficient for the channel in which the photon produces an 
electron in Landau level j and a positron in Landau level k, 
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and the sum is taken over all possible states for the electron- 
positron pair. Since pair production is symmetric with re- 
spect to the electron and the positron, R'jf. = R'f^j', for sim- 
plicity we hereafter use R'jf. to represent the combined prob- 
ability of creating the pair in either the state (jk) or (kj) 



12 



/ new 

■jk 



R 



+ R'kl )■ For a given channel {jk), the 



(i.e., R 

threshold condition for pair production is 

e > E,+Ek, 



(23) 



where e' — e sin ^ is the photon energy in the perpendicular 
frame and E„ = nieC^y^l + ifign is the minimum energy 
of an electron/positron in Landau level n (the energy of an 
electron/positron with the momentum along the magnetic 
field P|| = 0). In dimensionless form, the condition [Eq. (|23p ] 
can be written as 



2m„c^ 



i [Vl + 2/?Qj + x/l + 2/3Qfc] . (24) 



Note that Xjk satisfies 



Xll < X03 < 
X03 < Xll < 



(25) 



The first three attenuation coefficients (corresponding to the 
three lowest threshold levels a;oo, a:oi, 0:02) for both || and _L 
polarizations are given in Appendix[Bl Eqs. (fB6| - (|BT0| : see 
also .Dau ghcrtv & Harding (1983). Note that -Rx.oo — 0, and 
thus the first non-zero attenuation coefficient for _L polarized 
photons is actually -R^ oii ^i^it -Rx.oo- 

In our simulation a photon is typically created with 
X below the first threshold (a;oo or xq\, depending on the 
photon polarization). As long as x remains below the first 
threshold, R' — and the optical depth to pair production 
remains zero. As the photon propagates into the magneto- 
sphere and crosses the first threshold, R' > 0, Ar > 0, 
and r begins to grow. As it continues to travel outward, 
both T and the number of Landau levels available for pair 
production jmax and fcmax increase. Depending on the local 
magnetic field strength [Eq. (|lip ]. the photon may reach a 
large enough optical depth (r ~ 1) for pair production after 
crossing only a few thresholds (so that jmax and fcmax are 
small) or after crossing many thresholds (so that jmax and 
femax are very large). For "weak" magnetic fields (/3q <0.1) 
the optical depth increases slowly with Sph and it is valid to 
use the jmax,fcniax ^ 1 asymptotic attenuation coefficient 



for pair production (e.g.. lErbeilll96d ). 



R'n 



0.23 

/ 

ao 



, exp 



SxPqJ 



(26) 



which applies for both polarizations. For stronger fields, 
however, pairs are produced in low Landau levels, and the 
more accurate coefficients of Daugherty & Harding must be 
used. In Appendix lB2l we find that the critical magnetic field 
strength separating these two regimes is 



Be 



(27) 



[see Eq. (|B17|) ]. We also find that the boundary between the 
two regimes is very sharp: pairs are either created at the first 
few Landau levels (n < 2) for B > Bcrit or in very high Lan- 
dau levels for _B<_Bcrit, with very few electrons/positrons 
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Figure 5. Attenuation coefficients in the perpendicular frame 
(the frame where the photon is traveling perpendicular to the 
magnetic field), for both photon splitting, labeled by -L— >'|| -I- ||, 
and pair production, labeled by \\—^ e+e~ and -L— > e~^e~. The 
local magnetic field strength is B = Bq = 4.414 x 10^^ G. 



created in intermediate Landau levels. Therefore, in our sim- 
ulation we only consider the first three attenuation coeffi- 
cients for ||-polarized photons (iiy oo)^(| ou^Il 02) ^'^'^ the 
first two non-zero attenuation coefficients for _L-polarized 
photons (i?x,oi, ^1,02)- If the photon reaches the threshold 
for the (03) or (11) channel [whichever is reached first; see 
Eq. (|25|) ]. we use the asymptotic formula, Eq. (|26|) . The to- 
tal attenuation coefficient for pair production (as given by 
this approximation) is plotted in Fig. [5] for both || and _L 
polarizations at /3q = 1. 

We include photon spli tting in our simulations . Based 
on the kinetic selection rule (|Adleij|l97ll : IUsovll2002l . but see 
iBaring fc HardindlioOll ) . only the process -L— ^|| || is allowed. 
Therefore, for [[-polarized photons, the attenuation coeffi- 
cient for photon splitting is zero (i?,jl'P = 0). For _L-polarized 
photons we use the following formula, adapted from the nu- 
merical calculation of iBaring fc Harding! (|l997l ): 



[g{PQ,x) + Q.Qb] [0.25p(/3Q,a;) + 20] 



(28) 



where g{l3Q, x) = /3q exp(— O.Oa::'^). For x < 1, this expression 
reproduces the results of Baring & Harding to better than 
10% at both Pq < 0.5 and /3q 3> 1, while underestimat- 
ing the results at Pq = 1 by less than 30%. The -L— >[| [| 
attenuation coefficient for photon splitting is plotted in 
Fig. [5] at B = Bq. Because the attenuation coefficient R'^'^ 
drops rapidly with field strength for /Jq < 1, photon split- 
ting is unimportant for /Jq < 0.5 (e.g., IBaring fc Harding 
[200^. However, for _L-polarized photons propagating in su- 
perstrong fields Pq > 0.5, photon splitting is the dominant 
attenuation process: even though above the first threshold 
{x > xoi for _L photons) the attenuation coefficient for pho- 
ton splitting is much smaller than that for pair production, 
in superstrong fields the photon splits before reaching the 
first threshold (see Fig. [5]). 

In the simulation, whenever r > 1 or Tsp > 1 the pho- 
ton is destroyed (i.e., turned into a pair or two photons). 
More precisely, the photon should only be destroyed with 
probability 1 — exp(— r). But in practice we find that such 
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a refinement has a neghgible effect on the cascade result. If 
Tsp > 1 the photon splits into two. As a simplification we 
assume that each photon takes half of the energy of the par- 
ent photon fcf. iBaring fc Harding 1997i ): therefore, at the 
point of photon splitting a new photon is created with an 
energy 0.5eo and a weighting factor 2A.N^ (i.e, the simula- 
tion photon represents two actual photons). The new pho- 
ton is ||-polarized (for the -L— >|j || process) and is assumed 
to be traveling in the same direction as the parent pho- 
ton, Oph. If r > 1 the photon creates an electron-positron 
pair. For B < Bcrit ~ 3 X 10^^ G, the pairs are created in 
high Landau levels (see above), and we assume that the 
electron and positron each shares half of the photon en- 
ergy and travels in the same direction as the photon: thus 
^ym^c? = e/2 and the electron/positron's magnetic pitch 
angle is ^ = ij). This app roximatio n is v alid as long as 
(see iDaughertv fc Hardina' ll983l ). which accord- 
ing to ML07 is satisfied for _B<Bcrit. When _B>Bcrit, the 
electron and positron are created in low Landau levels (we 
choose the maximum allowed values, jmax,fcinax, since this 
channel dominates the total attenuation coefficient), with 
energies given by Eq. (|B5|I of Appendix IBll 

In the simulation we try to find the photon-magnetic 
field intersection angle at which pair creation occurs, i/'pair, 
to an error of less than 10%. If the error in i/'pair is too large, 
the electron and positron will be created in the wrong Lan- 
dau level and will emit too many or too few synchrotron 
photons (see Section [3^ . To accurately determine i/'pair we 
use the following procedure in our simulation: The photon's 
first full step. So, should be small enough that the probabil- 
ity of pair production is negligible at so but large enough 
that the probability grows rapidly with subsequent steps. 
At high fields a good choice for so is the location of the first 
non-zero threshold {x — xqo for || polarization or a; = xqi 
for _L polarization), since the attenuation coefficients are 
large enough to allow pair production in a distance much 
shorter than 1 cm. At low fields a good choice is the point 
where x/Bq — 1/20. At this point the mean free path for 
pair production is much larger than the gap height while for 
xPq = 1/10, e.g., the mean free path is much smaller than 
the gap height. Therefore, so is chosen such that it solves 

fa;oo , /3q > 1/20 and || polarization; 

xoi , /3q > 1/20 and _L polarization; (29) 

1/(20/3q), I3q< 1/20. 

Note that both x and Pq depend on so, so the value of so 
must be found numerically. Since simp ~ Sph/TZc (for small 
angles), this distance is approximately given by 

„2 



So 



:7^, 



2meC 



1 + 



20/?Q 



(30) 



In our simulation, for || polarizations (no photon splitting), 
the photon moves directly to so in one step. For _L polar- 
izations, the photon moves to so in 10 steps (with step sizes 
0.1 So), allowing for the possibility of photon splitting before 
reaching this point. In either case, once the photon reaches 
So it steps outward in the manner described at the beginning 
of this section. At high fields {B> Bait) we choose the step 
size to be Asph ~ 0.1so(xoi — 2;oo)/a;oo for || polarizations 
or 0.1so(a;o2 — a;oi)/xo2 for _L polarizations. At low fields 
(_B<i3crit) we choose Asph = O.lso. 

Sometimes the photon does not pair produce (or split) 



before exiting the magnetosphere. Conveniently, we do not 
have to track the photons out to the light cylinder to know 
whether pair production will occur. Once a photon reaches 
the xo3 or xn threshold, such that the asymptotic expres- 
sion for pair production Eq. (|26|l can be used (i.e., when 
B<Bciit), then the growth in optical depth depends "ex- 
ponentially" on xPq [since At cx exp{— 1/(2;/3q)}]. Because 
xPq oc Sph(ro,ph + Sph)~^^^ reaches a maximu m at Sph — 
0.4ro.p h and then rapidly decreases (cf. Hibschman fc Aronj 
l2001al ). we assume in our simulation that if the photon does 
not pair produce by 



Sph,n 



0.5ro,ph : 



(31) 



it will never pair produce and instead escapes the magneto- 
sphere. Here r-o,ph is the altitude of the photon at the emis- 
sion point. Note that Eq. (|31|l is also approximately valid 
for our treatment of non-dipole fields {TZc ~ R near the stel- 
lar surface), since once the photon has traveled a distance 
Sph — Sph, max it is iu the dipole regime (r > 27?). 



3.3 Propagation and photon emission of the 
secondary electrons and positrons 

3.3.1 Synchrotron radtatton 

In the corotating frame (the "lab" frame) the secondary elec- 
tron (or positron) is created with energy 'yrrieC? , pitch an- 
gle 'I' [with the corresponding Landau level n; see Eq. (|32p 
below] and weighting factor AA^^ (Section lS.lf) . For the pur- 
pose of tracking the synchrotron emission from the electron 
it is easier to work in the "circular" frame, the frame in 
which the electron has no momentum along the magnetic 
field direction and only moves transverse to the field in a 
circular motion. Note that this frame is in general different 
from the perpendicular frame (defined in Section [3]2} of the 
progenitor photon; only if the electron-positron pair is cre- 
ated exactly at threshold \x = Xjk', see Eq. (I24p ] are the two 
frames the same. The energy of the electron in the circular 
frame, E±_ = 'y^rUeC? , is related to that in the lab frame by 



71 



r 



+ cos2 = y/TTWon . 



(32) 



Note this expression also gives a relation between 7x and 
n; we shall use 7x and n interchangeably to refer to the 
electron's energy in the circular frame. 

In the circular frame i5x is radiated away through syn- 
chrotron emission on the timescale 



synch 



£x 



^synch 

5 X 10" 



7xmeC 



24 (-Y? 



— l)c^Wc 



(33) 



where Uc = eB jm^c is the electron cyclotron frequency and 
B\2 is the local magnetic field strength B in units of 10^^ G. 
This decay time is much shorter than other relevant cascade 
timescales (e.g., the timescale for B or TZc to change sig- 
nificantly, which is of order r/c >10~'' s, or the timescale 
for the emission of resonant ICS photons, discussed later 
in this section). Therefore, in the simulation the electron 
is assumed to lose all of its perpendicular momentum px 
"instantaneously" due to sync hrotron radiation, before mov - 
ing from its initial position (cf. IDaughertv fc Hardin3ll982l ). 
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The final energy of tfie electron once it reaches the ground 
Landau level (or p± — 0) is given by 



7|| 



(1-/3^ cos^*)"^/^ = 7/71 



(34) 



where /3 — \J\ — 1/7^ is the electron velocity. 

Since the synchrotron photon may carry an energy com- 
parable to E\_ of the parent electron, it is necessary to track 
the electron energy after each photon is emitted in order 
to obtain accurate synchrotron spectrum (this is in contrast 
to the case of curvature radiation discussed in Section 13.11 
where a large number of curvature photons can be emit- 
ted without significantly affecting the energy of the parent 
electron) . As a simplification, in the circular frame the syn- 
chrotron photons are assumed to be emitted isotropically in 
the plane of motion, such that no velocity kick is imparted 
to the electron; thus the frame corresponding to circular mo- 
tion of the electron does not change over the course of the 
synchrotron emission process. In other words, as the elec- 
tron loses its px_ , the Lorentz factors 7 and 7x decrease but 
7|l is constant, and Eq. (|34|l remains valid during the entire 
synchrotron emission process. 

We adopt the following procedure in our simulation: 
In the circular frame, the electron Lorentz factor 7^ drops 
from its initial value to 7^ = 1 in a series of steps; when 
7_L = 1 (i.e., n = 0) synchrotron emission stops. In each 
step one synchrotron photon is emitted, with an energy e± 
that depends strongly on the "current" value of y±. After 
the photon is emitted the energy of the electron is reduced 
by the amount A7X = e±/me(?. In the next step another 
photon is emitted with a new value of ex, and so on. 

In the simulation the photon energy e±_ of the syn- 
chrotron radiation is chosen in one of three ways, depending 
on the Landau level number n of the electron, (i) If the 
electron is created in a high Landau level (n > 3), the en- 
ergy of the photon is chosen randomly, but with a weight- 
ing based on the asymptotic synchrotron spect runfl (e.g., 
ISokolov fc Ternovlll968l : iHarding fc Preecelll98'iD 



dtde I 



where 



V3 



2tt 



fF 



+ 



ex 



G 



3 2t, 



(35) 



(36) 



is the characteristic energy of the synchrotron photons, 
f ~ 1 — ex/(7x77ieC'^) is the fraction of the electron's energy 
remaining after photon emission, F{x) — x Kf,/-i{t) dt, 

^ Thi s expression diff ers from the classical synchrotron spectrum 
(e.g., iRvbicki fc Ligh tman 1979.3 ™ ways: First, a factor 

of / = 1 — ex/(7x'TteC^) appears in several places throughout 
Eq. II35I I: when the photon energy is equal to the electron energy 
(ex = 7x™eC^ or / = 0) the asymptotic expression goes to zero. 
Second, a term containing the function G(x) appears in Eq. 1351 1. 
While such a term appears in the classical expressions for the 
radiation spectra of both Ij- and ^-polarized photons, in the clas- 
sical expression for the total radiation spectra these terms cancel 
out. However, when the quantum effect of the electron spin is 
considered there is an asymmetry between the perpendicular and 
parallel polarizations such that the G{x) term remains. 



and G{x) = xK2/-i{x) [cf. Eq. (ii) If n = 2, the 

energy of the photon is either that required to lower the 
electron to its ground state (n = 2 — > 0) or the first 
excited state (n = 2 — > 1), with a probability that de- 
pends on the local magnetic field strength. We do not use 
the exact transition rates for the n = 2 state here. In- 
stead, we use t he following simplified pre s cripti on, based 
on the results of Hcrold. Ruder. &: Wunnen (| 19821 ') (see also 
IHarding fc Preece 1987): If /3q < 1 the energy of the pho- 
ton is chosen to be that required to lower the electron to the 

first excited state, ex = meC^ ^^1 -I- 4/3q — ^1 -f 2j3Q^. If 

/3q > 1 the energy of the photon is randomly chosen to be 
that required to lower the electron to either the ground state 
[ex = meC^ (^^1 -I- 4/3q — 1^], 50% of the time, or the first 

excited state, 50% of the time, (iii) If n = 1, the energy of the 
photon is that required to lower the electron to its ground 
state, ex = rUeC? ^\/l -f 2/3q — 1^ . If the electron is not in 
the ground state after emission of the synchrotron photon 
(which could happen for the n = 2 and n > 3 cases discussed 
above, but not for the n = 1 case), 7x is recalculated and a 
new photon energy is chosen. 

The energy of the photon is transformed from the cir- 
cular frame into the "lab" frame using 



711 f 



(37) 



The photon carries with it the same weighting factor AA^'e as 
the secondary particle that emitted it. Because the photon is 
emitted in a random direction perpendicular to the magnetic 
field in the circular frame, in the lab frame the angle of 
emission (relative to the dipole axis) is approximately given 
by 



e 



ph 



X + * cos n , 



(38) 



where 11 is a random angle between and 27r, x is the angle 
between the local magnetic field and the dipole axis and is 
given by Eq. (fT6|l , and the pitch angle is given by Eqs. p2)l 
and 131: 



7i 



(39) 



For synchrotron radiation the polarization fraction is be- 
tween 50% and 100% polarized perpendicular to the mag- 
netic field (which is the exact opposite of th e curvature ra- 
diation case; see iRvbicki fc LightmanI [l979l ) . Therefore we 
randomly assign the photon a polarization in the ratio of 
one II to every seven _L photons (corresponding to a 75% 
perpendicular polarization). 



3.3.2 Resonant inverse Compton scattering 

Once the electron loses all of its perpendicular momentum, 
it moves along the magnetic field line in a stepwise fashion 
while upscattering surface thermal photons through RIGS. 
The step size As is related to AArics, the number of pho- 
tons scattered in each step, by 



As : 



cAArics 
dNmcs/dt 



(40) 
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In our simulation we choose AA''rics to be 

AJVrics = min (^1 , O.lR^^^^^^fl^^ . (41) 

In other words, AA^rics = 1 if the RICS process is efficient 
enough to produce at least one resonant photon within a 
distance of 0.17?; otherwise AA'rics is chosen so that the 
electron step size is As = O.IR. Using Eq. (|A6P from Ap- 
pendix [XJ we have 



As : 



AA^Rics 



7/3 



ag \ 1 



In- 



c/[7||(l-/3||)fcT] 



c/lT||{l-/3|| 



(42) 



where /^y — y^l — l/7y is the speed of the electron after it 
has completed synchrotron emission (so that p± = 0), and 
V'crit is the incidence angle with respect to the electron's 
trajectory of photons coming from the edge of the surface 
"hot spot" (see Fig. \KQ^ The mean energy of th e scattered 
photons is (;e.g.. lBeloborodov fc Thompsonll2007l ) 



(43) 



e = 7|| 1- 



v/1 + Wq 

and the energy loss of the electron in each step is given by 



A7||meC = -eAA/'rics ■ 



(44) 



In the lab frame the photon's angle of emission is ap- 
proximately given by 



9ph ~ X + — cos n 

7e 



(45) 



where 11 is a random angle between and 2n and x is the 
angle between the local magnetic field and the dipole axis 
[Eq. (|16| . Here, 7e is the final Lorentz factor of the elec- 
tron after emitting the photon; from Eq. (|43p. its value is 
approximately 



7e 



711 



\/l + 2/3q 



(46) 



While we include the l/7eCosn term in Eq. (|45() for com- 
pleteness, we find that it is not important for our simulation. 
This is true even at _B > Bq , where 7e is much smaller than 
the initial Lorentz factor of the electron and pair produc- 
tion c an occur almost immediate ly after the photon is scat- 
tered (jBeloborodov fc Thomps oni .2007 ). The extra distance 
traveled by the photons in order to pair produce when the 
photons are upscattered tangent to the local magnetic field 
(i.e., when 0ph = x is assumed) has a negligible effect on 
the overall cascade. 

In the superstrong field regime, the final polarization 
state of a photon upscattered th rough RICS is given by the 
results of iGonthier et~al] (|200(]| '). For B<Bq, both below 
and above resonance more _L photons are produced than || 
photons, at a ratio of ~ 3 : 1. The same situation occurs 
for B > Bq below resonance; above resonance, however, the 
situation reverses and more || photons are produced than _L 
photons. We therefore assign the photons a polarization in 
the ratio of one || to every three _L photons for B < Bq, and 
a polarization in the ratio of one || to every _L photon for B > 
Bq (based on the assumption that approximately 50% of the 
photons are slightly below resonance and 50% are slightly 
above). In practice, however, we find that the cascade does 



not depend sensitively on the initial photon polarization. At 
low fields (-B < 3 X 10^^ G) the polarization has no effect on 
the cascade, since the asymptotic attenuation coefficient for 
pair production is used; at high fields a _L photon is split into 
two II photons before it can pair produce, and the resulting 
cascade is not much different from the cascade of a single | 
photon with twice the energy. 

In our simulation we consider thermal photon emission 
from three types of surface "hot spots" (see lPons et al.ll2007l 
for a review of neutron star surface temperatures in strong 
magnetic fields): a large cool spot, = 0.3 and 6'apot = 7r/2, 
representing emission from the entire surface of a neutron 
star; a mid-sized warm spot, Te = 1-0 and Sapot = 0.3; and 
a small hot spot, Te = 3.0 and Ospot = 0.1, representing 
emission from a heated polar region 1 km across (which, 
though small, is still significantly larger than the polar cap 
region, unless P < 0.01 s). We also consider the case where 
ICS has no effect on the cascade, which we find occurs for 
Te < 0.1 (neutron stars too cold), Sspot ^ 0.01 (hot spots too 
small), or ro > -R(l -f- 2^spot) (particles injected too far away 
from the surface; this is most relevant for photon-initiated 
cascades discussed in Section [3l4l below) . 



3.4 Cascades initiated by a primary photon 

In the second version of the simulation, a photon is created 
with energy eo at the position (ro.ph, ^o,ph). We typically 
choose ro,ph = R and 6*0, ph ~ dcap (cf. Section [5TT]l . since the 
resonant ICS photon density is largest at r ~ _R; however, we 
are also interested in photons emitted at a higher altitude 
(e.g., for surface field strengths Bp > 10^'' G the behavior of 
the cascade with ro,ph = R and with ro,ph ~ 3R are very 
different; see Section We set AA^'e = 1, such that each 
photon in the simulation represents exactly one photon in 
reality; we can later multiply the simulation results by A^o of 
Eq. @ if we wish to compare cascades dominated by RICS 
and by curvature radiation (see Section [2. 2p . The photon is 
injected tangent to the magnetic field [Qph = Xi cf. Eq. (|17|l ]. 
since we find almost no difference in the final photon or pair 



spectra if we add a beaming angle I/7 



10" 



-lO-". As 



was discussed in Section [3.3.21 resonant ICS photons have a 
polarization ratio || to _L of approximately 1:3 for B < Bq 
and approximately 1:1 for B > Bq. In Section [4.11 results 
for the photon-initiated cascades, we choose the initial pho- 
ton to be polarized perpendicular to the magnetic field, to 
create a cascade with particle multiplicities as large as pos- 
sible. Our results therefore represent upper bounds on the 
actual cascade multiplicities. The actual cascade should not 
differ greatly from that presented in Section [4. II however, as 
cascades initiated by photons polarized parallel to the mag- 
netic field are only slightly lower in particle multiplicity and 
are qualitatively similar in spectral shape. Once the initial 
parameters of the photon have been chosen, the simulation 
proceeds in the exact same way as described in Sections l3.2f - 
13.31 the photon steps outward in a straight line from the 
point of emission until its optical depth is large enough to 
pair produce or split, etc. 

Note that due to the discrete, random nature of the 
synchrotron emission and the small number of particles in- 
volved in the cascade, photon-initiated cascades will have 
photon and pair spectra that are coarse and that vary be- 
tween simulation runs. In order to smooth/average the spec- 
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tra to some extent, we modify the synchrotron emission pro- 
cedure of Section |3]3] for secondary particles in high Landau 
levels (n > 3). In every step 10 photons are emitted, each 
with a weighting factor of 0.1 AN^^ (rather than one photon 
with a weighting factor of AN^, as before). Each photon 
has a different energy e±^i [chosen randomly according to 
Eq. (|35|) ]. so that the total energy lost by the secondary 
particle becomes 7j_meC'^ = 0.1 X^I^i ^-L,*- We do not apply 
this procedure to the synchrotron emission from secondary 
particles in Landau levels n = 1 or n = 2, as it would not 
gain anything; each of the 10 photons emitted would have 
the same value of ex i- 



4 RESULTS 

In this section we present the results of our simulations of 
photon- and electron-initiated cascades (Sections 14.11 and 
14.21 respectively), for a variety of different surface field 
strengths, rotation periods, field geometries, and initial en- 
ergies of the primary particle. For each type of cascades, we 
present the "final" spectra of the cascade photons and pairs 
as they cross the light cylinder and escape from the magneto- 
sphere. For the electron-initiated cascades we also show the 
spectra at several intermediate stages (i.e., the spectra of all 
photons and pairs that cross the height r = 1.2/?, 2R, 5R, 
etc.). The photon spectra are plotted over the energy range 
10 keV-1 TeV, since for energies < 1 keV the thermal pho- 
tons dominate the spectra while above ~ 1 TeV fewer than 
one photon is produced per primary electron. We are par- 
ticularly interested in the pair multiplicities, i.e., the total 
number of cascade electrons -I- positrons produced per pri- 
mary particle. We use n_E to denote the number of electrons 
and positrons per "primary" photon and Ne to denote the 
number per primary electron; the two multiplicities are re- 
lated by 



Ne = No X ue , 



(47) 



where A^o is the number of photons produced by the primary 
electron (see Section [2. 2 |l . From our numerical results we in- 
fer various empirical relations for each cascade; quantitative 
arguments for the validity of several of these relations are 
given in Appendix [Cl 

We first present our results for photon-initiated cas- 
cades (see Section l3.4p . as they are simpler and aid us in 
our discussion of the results for the full cascade (initiated 
by a primary electron). 

4.1 Results: photon-initiated cascades 

Our results for photon-initiated cascades are presented in 
Figs. ISHSl We consider primary photons with energies in 
the range of 10^-10^ MeV; for Bp, 12 = 1-1000, the primary 
electron should emit very few photons (via either resonant 
ICS or curvature radiation) above this energy range (see 
Section [2} . Unless otherwise stated, the primary photon is 
emitted from near the surface, in the direction tangent to 
the last open field line. Thus the radius of curvature near 
the point of emission is TZc — 9 x IO'^Pq ''^ cm for dipole fields 
[Eq. (HI]. 

We find significant differences in the behavior of the 
cascades at magnetic field strengths below and above 



Bcrit ~ 3 X 10^^ G [Eq. (|BT7)| ]. At low fields B<Bcrit, 
the primary photon can pa ir produce if [Eq. (|C2|) : see also 
iHibschman fc Aronsir2001al | 



eo > En 



3000B~12 7^8 MeV, 



(48) 



where TZs is the radius of curvature TZc in units of 10 cm, 
evaluated at the surface along the last open field line. Strong 
cascades, where more than one electron-positron pair is pro- 
duced, typically occur at energies ~ 10 times emin. For eo in 
the range from emin to ~ lO'' MeV, we find that the mul- 
tiplicities of photons and e^e~ particles produced in the 
cascade are 



eo 



500 MeV 



7^8 



and 



riE 



(49) 



(50) 



10^ MeV 

respectively. These results are (largely) independent of the 
hot spot model used. When ICS is inactive, the cascade elec- 
tron/positron has final energy (after it has finished radiating 
synchrotron photons) extending from [Eq. (|C3|) ] 



En 



0.1-Bp,i2eo 



(51) 



2En 



0.1i3p,i2eniinn£; In 



(for the first pair produced) down to ~ 0.1-Bp,i2emin for 
the lowest-energy pairs, and the total energy of the pairs 
is [Eq. ^] 

'0.0_75eo\ (32) 

mill / 

When ICS is active from a hot spot (Section l3.3|l . the num- 
ber of pairs produced does not change, since the photons 
produced through ICS at these field strengths have energies 
~ Bp^i2Tff^ MeV [Eq. ((H])] and can not pair produce. The 
total pair energy Stot decreases, however, since the ICS pro- 
cess transfers energy from the pairs to photons. Although 
resonant ICS is most important for electrons and positrons 
at 7crit — (.c/kT [Eq. ((7]); see Section [2.2j . we find in these 
cascades that all electrons and positrons with energies in the 
range of 



-Erics ~ (0.3 - 30)7critmeC^ 

~ (20 - 2000)Bp,i2r6"^ MeV 



(53) 



are strongly affected. Thus hot surface spots with higher T 
tend to lower £tot more. As expected, we find that photon 
splitting does not affect the cascade at these field strengths 
(see Section 13. 2p . The photon and pair cascade spectra for 
Bp, 12 = 1 are shown in Fig. [6] both when ICS is inactive and 



when ICS is active from a "warm spot" (Te 



spot 



= 0.3). 



At high fields (_B>_Bcrit), a primary photon injected 
from the surface will pair produce when 



eo > e,, 



200 7^8 MeV, 



(54) 



largely independent of field strength. When ICS is inactive, 
almost all of the cascade energy resides in the pairs; i.e., 
Etot — eo. The pair cascade will be very weak regardless of 
photon energy, with ue < 10 and tIj = or 1 (i.e., at most 
one photon escapes the magnetosphere without pair produc- 
tion). This is because the pairs are produced exclusively 
through the (jk) = (00) or (01) channel (see Section[32|, so 
that at most one synchrotron photon is emitted per pair. For 
Bp, 12 ^ 20, photon splitting causes all pairs to be produced 
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Figure 6. The final photon and pair spectra of photon-initiated cascades for surface magnetic fields Bp, 12 = 1- The NS spin period is 
Po = 1 and a dipole field geometry is adopted. In the upper panels, ICS is assumed to be inactive, while in the lower panels, ICS from a 
hot spot with Tg = 1, (9spot = 0.3 is included in the simulation. The primary photon is injected from the surface and has an energy of 
10^ MeV (left panels) or 10^ MeV (right panels); for photons with energy 10"^ MeV, no cascade is initiated. The spike in the pair spectra 
of each panel represents the electron-positron pair produced by the primary photon. The spectra in the top panels (where ICS is inactive) 
are nearly identical to the spectra generated, e.g., by a photon injected at ro,ph = 3-R above a star with surface field Bp, 12 = 3^ = 27, 
such that the local field strength at the injection point is B = 10^^ G [Eq. mil l]. 



with (jk) = (00), such that the cascades are even weaker: 
n_B < 4 and = 0. "When ICS is active, both and he can 
be larger, though not as large as would be predicted by an 
extrapolation of Eqs. (|49l) and (|50[) to high fields. In order 
for ICS to aflect the cascade, however, the primary photon 
must have an energy 

eo>70Bp,i2Tg-' MeV; (55) 

the energy of the electron/positron produced by this photon, 
^'max ~ (0.1-0.5)eo [Eq- (|C9[) ]. must be larger than the min- 
imum energy at which ICS is effective, ~ 0.37crit'7ieC^ (see 
Appendix lCl|) . Note that this energy is approximately equal 



to Erics , the energy of a typical ICS photon upscattered by 
the primary electron at high fields [Eq. ((SJ]; therefore, at 
high fields a typical ICS photon is able to initiate a weak 
cascade. We find that for Bp < 0.5Bq a significant fraction 
(~ 20%-60%) of the total cascade energy eo resides in the 
photons; as in the low field case [Eq. ((52])], this fraction de- 
creases as either eo or -Bp, 12 increases. For Bp>0.5BQ, even 
ignoring photon splitting, the total photon energy fraction 
is very low, < 10%; but with photon splitting included, it is 
almost negligible, < 1%. The photon and pair cascade spec- 
tra for -Bp, 12 = 10 and 100 when ICS is active from a "warm 
spot" (Te = 1, e^pot = 0.3) are shown in Fig.0 
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Figure 7. The final photon and pair spectra of photon-initiated cascades with active ICS, for surface magnetic fields Bp, 12 = 10 
(upper panels) and 100 (lower panels). The pulsar spin period is Pq = 1, a dipole field geometry is adopted, photon splitting -L— ^|| || 
is active, and thermal photons are emitted from a hot surface spot with Ta = 1, 6spot = 0.3. The primary photon has an energy of 
10^ MeV (left column of panels), 10* MeV (middle column), or 10^ MeV (right column). The spikes in the pair spectra of several panels 
represent the electron-positron pair produced by the primary photon. Note that there are actually two such pairs in the bottom left panel 
(Sp,l2 = 100,60 = 10'^ MeV), as the primary photon has split in that case; while in the upper right panel (Bp, 12 = 10, eo = 10^ MeV) 
the electron and positron have significantly different energies from each other and so are represented by two shorter spikes. Also note 
that photon spectra do not appear in the lower panels. This is due to a combination of weak synchrotron emission and efficient pair 
production near the surface of a Bp, 12 > 100 neutron star; very few secondary photons are created, and none of them survive to escape 
the magnetosphere. 



As discussed in Section [3.31 we consider three hot sur- 
face spot models for active ICS: a "cool" Te = 0.3, Sgpot = 
7r/2 spot; a "warm" Te = 1, Sspot = 0.3 spot; and "hot" 
Te = 3, Sspot = 0.1 spot. At low fields we find that the 
only effect the various hot spot models have is to lower £tot 
relative to the total cascade energy eo (see above). At high 
fields, the cascades due to warm and hot spots are similar in 
multiplicities and ue for energies eo < IC* MeV; but for 
energies eo > 10^ MeV, the hot spots give multiplicities ~ 3 



times larger than warm spots; see Fig. [S] [from Eq. (|55p . at 
Bp,i2 ~ 1000 a strong cascade requires Te > 3]. Cool spots 
give much smaller multiplicities (factors of > 10 smaller) 
than warm or hot spots regardless of the primary photon 
energy. 

At low fields, the multiplicities of photons and e'^e~ 
particles produced per primary photon, and ue, depend 
on TZc to the ( — 1) power [Eqs. (|49)) and ((50)) ]. such that 
cascades in non-dipole magnetospheres can be much larger 
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Figure 8. The final photon and pair spectra of photon-initiated cascades with active ICS, for surface magnetic fields Bp, 12 = 10 (left 
panel), 100 (middle panel), and 1000 (right panel). The NS period is Pq = 1, a dipole field geometry is adopted, and thermal photons 
are emitted from a hot spot with Tg = 3, Sspot = 0.1. The primary photon has cq = 10^ MeV. The spikes in the pair spectra of the left 
and right panels represent, from shortest to tallest, one, two, or four electrons/positrons. 



than in dipole magnetospheres. At high fields, we find that 
this dependence on curvature radius is much weaker: ue oc 
TZc^^"^ at Bp ~ 0.5-Bq and ue oc TZZ^^'^ power at larger 
fields. Thus for high fields, the cascades in non-dipole and 
dipole magnetospheres are of similar sizes. The photon and 
pair cascade spectra for -Bp, 12 = 1, 10, and 100 and a TZc = R 
non-dipole magnetosphere are shown in Fig. (9] 

Two aspects of the cascade depend strongly on the al- 
titude (»'o,ph) at which the primary photon is injected: the 
local magnetic field strength, B oc r~'^, and the effectiveness 
of ICS, which is completely negligible for ro,ph > 3i? (regard- 
less of the temperature and size of the hot spot). The other 
cascade parameters have a much weaker dependence on al- 
titude (e.g., radius of curvature TZc oc y^ro^). We find that 
the photon and pair spectra for primary photons injected 
at ro.ph > R are very similar to the spectra for photons 
injected at the surface, as long as ICS is inactive and the 
local magnetic field strengths are the same in both cases. 
For example, the spectra for r-o,ph = 3R and Bp, 12 ~ 27 
[such that the local field strength at the point of injection is 
B12 = 1; Eq. (|lip ] is nearly identical to the spectra given in 
the first row of Fig. [6l where ro.ph ~ R, Bp, 12 = 1, and ICS 
is inactive. 

We can use the pair multiplicity per primary photon, 
Ue (eo), obtained from our numerical simulations to estimate 
the pair multiplicity per primary electron, Ne = NqUe, 
when ICS is the dominant cascade emission process. For 
eo and A'o, we use the expressions for the typical energy 
ERICS and total number ~ 10B~l2Po^ ^tI''^ of resonant ICS 
photons upscattered by the primary electron [Eqs. ((8| and 
((91; see Section [2.2) . The results are shown in Table [1] Note 
that although erics is independent of the acceleration model 
used, the number of upscattered photons, A^o, as given by 
Eq. ((9| is applicable for only inner gap accelerators with 



space-charge-limited fiow; inner vacuum gap accelerators, 
for example, would yield A'o about 20-100 times lower. Also 
note that, for a given acceleration model, an accurate de- 
termination of Ne requires that rather than just setting 
eo = tRics, a distribution of energies be used which takes 
into account resonant (and possibly non-resonant) scattering 
away from the thermal peak (i.e., 7 7^ 7crit): although fewer 
photons are upscattered at energies greater than erics , these 
photons can have an important effect on the total multiplic- 
ity of pairs produced, since ue grows approximately linearly 
with photon energy eo; the eo > eRics photons are especially 
important for cascades where ue (eRics) = (e.g., cascades 
with Bp, 12 = 10 and Po = 1). 

4.2 Results: electron-initiated cascades 

Our results for electron-initiated cascades are presented in 
Table [2] and Figs. I10H15I In our simulation, the primary 
electron is emitted from the surface along the last open field 
line {Oo = Scap). We consider the cases of 70 = 2 x 10^ and 
4 X 10^ in dipole magnetospheres, as well as the case of 70 = 
2 X 10^ in non-dipole TZc = R magnetospheres, as discussed 
in Section 12.11 Although a larger initial primary energy 70 
gives rise to more cascade particles with a larger total energy, 
the behavior of the cascade at 70 > 4 x 10^ is qualitatively 
similar to that at 70 = (2-4) x 10^ (or 70 = 2 x 10*^ and 
TZc = R) and key quantities such as cascade multiplicities 
and energies can be extrapolated from our results. For most 
of this section we "turn off" ICS in our simulation and allow 
photon splitting only through the -L— || mode; at the end 
of this section we discuss how changing these simulation 
parameters affects the cascade. 

In Table [2] we list some key quantitative results of our 
simulations: For each cascade (characterized by the spin pe- 
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Figure 9. The final photon and pair spectra of photon-initiated cascades for a non-dipole magnetosphere with local radius of curvature 
TZc = ^i, for surface magnetic fields Bp 12 = 1 (left panel), 10 (middle panel), and 100 (right panel). Here, the primary photon has 
eo = 10** MeV and ICS is active from a hot spot with Tg = 1, 0spot = 0.3. 



Table 1. Pair multiplicity when resonant ICS is the dominant photon emission process of the primary electron. For simplicity, all 
photons are assumed to be upscattered from the thermal peak (7crit; see text). Each entry in the table gives Ne = Nque for a different 
surface magnetic field strength Bp, radius of curvature (either Tic = R or a dipole field curvature with the pulsar spin period Pq is 
specified), and hot spot temperature T and size ffspot- The pair multiplicity is zero for Bp, 12 = 1 and the hot spot models used here, 
even when TZc = R; we therefore omitted these entries from the table. 
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riod for dipole fields or the curvature radius for multipole 
fields, the surface field strength Bp and the primary electron 
energy 'yorUeC?), we give (the final energy of the pri- 

mary electron when it escapes the light cylinder), etot (the 
total energy of the cascade photons), ftot (the total energy of 
the secondary e^e~ pairs), and Ne (the multiplicity of e'''e~ 
pairs produced per primary electron). Note that the total 
cascade energy must satisfy 7omeC^ = -|- £tot + etot- 

We find that Ne is largest for cascades with strong sur- 
face fields, short rotation periods, multipole geometries, or 
large initial energies for the primary electron, but that re- 
gardless of cascade parameters the particle multiplicity sat- 
urates at Ne ~ 10*. Increasing Bp or 70, or decreasing P or 
TZc, tends to increase the ratio /b = £tot/{£tot + ftot); i.e., 
under these conditions a larger fraction of the "secondary" 
energy, 7omeC^ — 7/meC^, is transfered from the photons to 
the pairs. At low fields only a small fraction of the secondary 
energy is held by the secondary pairs (e.g., for -Bp, 12 = 1, 
/b <0.05), but at high fields this fraction is typically above 
50%'(e.g„ for Bp,i2 = 1000, /b>0.8). The average energy 



of a secondary electron/positron, E = £tot /Ne, is directly 
proportional to the total cascade energy (i.e., E oc 70), but 
depends only weakly on Bp, P, and TZc (e.g., E is approx- 
imately the same for TZc ~ R and 70 = 2 x 10^ as in the 
dipole case for 70 = 2 x lO'^). 

Figure [TD] shows the secondary pair multiplicities Ne as 
a function of 70, for various field strengths and periods/field 
geometries. We find that for dipole geometries Ne is a strong 
function of both 70 and the ratio BpP~^ (i.e., the polar cap 
voltage), but depends very weakly on either P or Bp alone. 
For example, the Ne versus 70 curves for -Bp, 12 = 1, Po = 1 
and -Bp, 12 = lOO, Po = 10 are nearly the same, as is shown in 
Fig. 1101 Assuming that a dense secondary pair plasma is a 
necessary ingredient in pulsed radio emission from neutron 
stars (see Section [T|, we can use our pair multiplicity results 
to estimate the conditions for pulsar "death" — the condi- 
tions under which a NS is no longer active as a radio pulsar: 
For a particular set of cascade parameters Bp, P, and TZc, 
we define 7death to be the value of 70 at which on average 
only one electron/positron is created per primary electron; 
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Table 2. Energies and multiplicities for a cascade initiated by a single electron. Listed are the initial and final energy of the primary 
electron, -yorrieC^ and "ffmeC^, respectively; the total energy of the cascade photons, etot; and the multiplicity and total and average 
energies of the secondary electrons and positrons, A'^;, £tot, E = £tot/NE, respectively. Each cascade is specified by the magnetic field 
strength, field geometry (curvature radius TZc or spin period in the case of dipole fields), and cascade energy (7omeC'^). 
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i.e., 

Ne (70 = Tdoath) = 1 • (56) 

The value of 7dcath changes very little if we alter the critical 
value of Ne in Eq. (|56|l by a factor of ~ 10, because of 
the steep dependence of Ne on 70 in this region. Therefore, 
although it is unknown exactly what value Ne must have 
for a pulsar to be active, 70 = 7dcath is a good predictor of 
pulsar death regardless. Using Fig. [10] we find empirically 
that 

7death ^ 1.5 X 10'B-l'/'7^^/' . (57) 

1/2 

For dipole fields, 7^8 = 0.9P(, ' , and we can write 
7doath =i 1.4 X WB-\!^''pI''' 

( 'I? death \ ^''^ ( e4'dcath\ 

where $cap is given by Eq. ^ and 

$dcath = 7 X 10^" V. (59) 
Therefore, for dipole fields, the death line is given approxi- 



mately by $cap = <&dcath, or 

^'o = B% ■ (60) 

for -Bp, 12 = 1-1000 this is nearly the same as the death line 
depicted in Fig.[T](Po = 0.65^''^^). 

An example of how the cascades develop spatially is pre- 
sented in Fig. 111! which shows the number of photons and 
secondary pairs as a function of the radius at which each par- 
ticle is created, for surface fields strengths -Bp, 12 = 1, 44.14 
(i.e., -Bp = -Bq), and 1000. The right two panels show the 
difference in the behavior of the cascade when the local mag- 
netic field strength B is above or below -Bcrit — 3 x lO'^^ G 
[see Eq. (I27|l ]. i.e., below or above r ~ 2.5-R for -Bp = Bq 
and r ~ 7i? for Bp,i2 = 1000 [see Eq. For B> Bcrit, 

pair production by curvature photons is very efficient; for 
the conditions depicted in Fig. 1111 (Pn = 1 and 70 = 2 x lO'^) 
approximately one pair is created for every two photons 
emitted by the primary electron. This is because about half 
of all the photons created at low altitudes lie above the limit 
emin ~ 200 MeV of Eq. ^ (see Fig. [TUl. Photon splitting 
has very little effect on the curvature photons, since the ma- 
jority (7/8) of these photons are [[-polarized. There is no 
synchrotron radiation for B > 0.5-Bq, since the [| photons 
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Figure 10. Secondary electron/positron multiplicities Ne as a function of the initial Lorentz factor of the primary electron 70, for 
Pq = 1 a-nd a dipole geometry at several different surface field strengths (left panel), and for Sp,i2 = 100 at several different periods/field 
geometries (right panel). 



will pair produce through the (jk) = (00) channel and _L 
photons will split into two || photons before pair producing 
(see Section[3]2} . For B < 0.5Bq, pair production dominates 
over photon splitting, such that a few synchrotron photons 
are emitted from electrons and positrons in the n = 1 Lan- 
dau level. For B < -Bcrit, electrons and positrons are created 
in higher Landau levels, such that many synchrotron pho- 
tons are produced by each electron or positron (~ 10; see 
Appendix EH. The minimum photon energy for pair pro- 
duction emin oc B~^ [Eq. (|48|) ] grows with radius while the 
energy of the typical curvature photon ecR oc 7^ [Eq. (|10p ] 
and synchrotron photon esR oc ecR (Appendix [Cl) fall with 
radius; therefore the number of electron-positron pairs cre- 
ated per photon decreases rapidly with radius. 

Figures and [T51 show the final spectra (i.e., the spec- 
tra as measured at the light cylinder) of photons and pairs, 
as well as the spectra of curvature photons emitted by the 
primary electron, for a variety of magnetic field strengths, 
spin periods, and cascade energies. Figure [14] shows cumu- 
lative photon and pair spectra at various magnetosphere 
radii. These spectra are generated by recording the energy 
of any photon, electron, or positron which reaches various 
radii (such as r/R = 1.2, 2, 5, 20). 

Figures [T^14l show many of the same trends as are seen 
in Tableland Fig. 1111 For example, the multiplicity Ne in- 
creases with increasing Bp, 70, , or TZc^ ■ But there are 
also several new features. First, increasing Bp or 70, or de- 
creasing P or TZc, tends to increase the maximum energy 
and decrease the minimum energy of the cascade pairs, so 
that the pair spectrum becomes broader in energy (see Sec- 
tion [C2]). Second, the cumulative photon and pair spectra 
at a given radius r (Fig. I14p are nearly independent of Bp, 
as long as the local field strength at that radius is larger 
than -Bcrit (e.g., the spectra at r = 2R are nearly the same 
for Bp, 12 = 44.14 and 1000). Third, under certain condi- 
tions, synchrotron radiation dominates the cumulative/final 
photon spectrum at low energies e ^ 1 MeV (e.g., in the top 
middle and top right panels of Fig. 113^ : For a given alti- 



tude r, the ratio of the number of synchrotron photons to 
the number of curvature photons at low energies is largest 
for large 70, small P, or small TZc- The ratio of synchrotron 
photons to curvature photons is also largest for Bp close to 
Bcrit, such that the point of maximum synchrotron radia- 
tion, B ~ Bcrit (see Fig. Illf) . occurs at a low altitude. Note 
that for e < 0.1 MeV the cumulative and final photon spec- 
tra all have power-law shapes with index F = 2/3 (where 
dN/de = (~^), regardless of whether synchrotron radiation 
from the secondary pairs or curvature radiation from the 
primary electron dominates; this is because both processes 
have spectra that depend on F{x)/e oc e~'^/^ at low energies 
[i.e., for x < 1; see Eqs. (fTS)) and pS^ ]. 

Figure [TS] shows the effect of resonant inverse Comp- 
ton scattering on the final photon and pair spectra. Both a 
"warm" hot spot (Ts = 1, 6'spot = 0.3), and a "cool" spot 
(Te = 0.3, Sapot = ■'r/2) are considered, as well as the case 
where ICS is inactive. For the cascade parameters adopted 
in our simulations, the results for "hot" hot spots (Tq = 3, 
^spot = 0.1) are nearly the same as for warm spots, and so 
are not shown. This is because the increased cascade activ- 
ity due to the larger T almost exactly cancels the decreased 
activity due to the smaller maximum photon-electron inter- 
section angle ■f/'crit). We find due to resonant ICS, a major- 
ity of the electrons and positrons with energies in the range 
.Erics ~ 20-2000Bp,i2rjri MeV [Eq. JSnj] radiatively cool 
to below that range. This has a strong effect on the pair 
spectra if the average electron/positron energy E lies in this 
range (which occurs, e.g., for Bp,i2< 44.14 when Po > 0.1 
and 70 < 4; see Table [2]). In general, we find that including 
RICS tends to make the cascade pair energy distribution 
narrower. On the other hand, ICS has only a minor effect 
on the photon spectra; its effect is moderate when E ~ 20- 
2000Bp,i2rg"^ MeV and Etot > etot (which only occurs when 
TZc = R and Bp_i2 = 10; see Table[5}. Regardless of the cas- 
cade parameters, resonant ICS by the secondary particles 
has very little effect on the multiplicity of photons or 
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Figure 11. The number of photons and electrons + positrons as a function of the radius where they are created, ro,ph or ''Oi for 
Bp,i2 = 1 (left panel), 44.14 (middle panel), and 1000 (right panel). The neutron star spin period is assumed to be Pq = 1 and a dipole 
field geometry is used. The -L— >|| || photon splitting mode is active, and the primary electron has 70 = 2 X lO'^; the effects of ICS are 
not included. The curve labeled "e+e~ pairs" shows where the secondary electrons and positrons are created, "Photons" shows where 
the photons that escape the magnetosphere (i.e., that do not split or pair produce) are created, "Curvature" shows where the curvature 
photons are emitted by the primary electron (which continues in a similar manner beyond the graph out to the light cylinder), and 
"Synchrotron" shows where the synchrotron photons are emitted by the secondary pairs. 
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Figure 12. The final photon and pair spectra of electron-initiated cascades. The primary electron has Lorentz factor 70 = 2 X 10^ (the 
left two panels) and 4 X lO'' (right two panels), and the surface magnetic field strengths are -Bp, 12 = 1, 44.14, and 1000. The neutron 
star spin period is Pq = 1 and a dipole field geometry is adopted. The secondary photon spectra are shown in the top two panels; the 
pair spectra are shown in the bottom two panels. The curve labeled "Curvature photons" in each of the top two panels (the dot-dashed 
line) shows the spectrum of curvature radiation emitted by the primary electron, which is the same for all field strengths. 
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Figure 13. The final photon and pair spectra of electron-initiated cascades for dipole magnetospheres with Pq = 10 (left panels) and 0.1 
(middle panels), and non-dipole magnetospheres with TZc = R (right panels). The surface field strengths are i3p,i2 = 1, 44.14, and 1000. 
The Lorentz factor of the primary electron is 70 = 2 X lO'^ for the dipole magnetospheres (left and middle columns) and 70 = 2 X 10^ for 
the non-dipole magnetosphere (right column; see Section [2J. The photon spectra are shown in the top three panels; the pair spectra are 
shown in the bottom three panels. The curve labeled "Curvature photons" in each of the top three panels (the dot-dashed line) shows 
the spectrum of curvature radiation emitted by the primary electron, which is the same for all field strengths. Note that the top right 
panel, corresponding to the photon spectra when TZc = R, has a different vertical (and horizontal) scale than the top left and top middle 
panels. 



pairs Ne created in the cascade, so it does not alter the 
pulsar death line given by Eq. (|58p . 

Figure [16] shows the effect of photon splitting on the 
final photon and pair spectra. We consider both the case 
where photon splitting is allowed only in the -L— >|| || mode 
(as discussed in Section 13.21 and in agreement with the se- 
lection rule of Adl er 1971; Usov 2002), and the case where 
both _L and || photons are allo wed to split (as is sug- 
gested in lBaring fc Harding||200ll ). For the cascade param- 
eters adopted in Fig. 1161 we find that when photon split- 
ting is "turned off" completely, the spectra are nearly in- 
distinguishable from the case where photon splitting is al- 
lowed only in the -L— || mode. In general, the effect of pho- 
ton splitting when both polarizations are allowed to split is 
strongest for large Bp and at high energies (e.g., large 70), 
such that the quantity x/Sq is large [see Eq. H28[) ]. At al- 
titudes where the local field strength is B > 0.5Bq (which 
occurs, e.g., at r <4i? for Bp, 12 = 1000), most photons con- 
tinue to split until they reach emin [Eq. H48p or H54p ] , and very 
few pairs are created. However, once the cascade reaches an 
altitude such that B < 0.5Bq, photon splitting has very 



little effect on the cascade. Because a majority of pairs are 
produced at altitudes above B ~ 0.5Bq (see Figs. [TT] and 
I14p . the pair spectrum is not strongly affected by photon 
splitting even when both polarizations are allowed to split 
and Bp,i2 = 1000. 



5 DISCUSSIONS 

We have presented numerical simulations of pair cascades in 
the open-field line regions of neutron star magnetospheres, 
for surface magnetic fields ranging from Bp = 10^^ G to 
10^^ G, rotation periods P = 0.1-10 s, and dipole and 
more co mplex field geometries. Compared to previous stud- 
ies (e.g.. iDaughertv Hardind Il982l: ISturner et all Il995 



Daughertv fc Hardind Il996l : jHibschman fc Aronj l2001b : 
Arendt fc Eilekll2002ri , which were restricted to weaker mag- 
netic fields (B < a few x 10^ G) and dipole geometry, we 
have incorporated in our simulations additional physical pro- 
cesses that are potentially important (especially in the high 
field regime) but were either neglected or crudely treated 
before, including photon splitting with the correct selection 
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Figure 14. The cumulative plioton and pair spectra at magnetosphere radii r = 1.05-R, 1.2R, 2R, 5R, 20R, and rLC (denoted "Light 
cylinder"), for Bp, 12 = 1 (left two panels), 44.14 (middle panels), and 1000 (right panels). The NS spin period is Pq = 1 and a dipole 
field geometry is used, and 70 = 2 X lO'^. The photon spectra are shown in the top three panels; the pair spectra are shown in the bottom 
three panels. For a given magnetosphere radius r, every cascade particle that is created below r and survives until reaching r is recorded 
in the spectrum. 




Figure 15. The ofi'ect of resonant inverse Compton scattering on the final photon and pair spectra. Two surface hot spot models are 
considered: a "warm" spot (Tg = 1, Sspot = 0.3) and "cool" spot (Ts = 0.3, fspot = '''/2). The case where ICS is inactive is also shown. 
The other cascade parameters are; Bp, 12 = 44.14, Pq = 1 (a dipole field geometry is used), and 70 = 4 X 10^. The photon spectra are 
shown in the left panel and the pair spectra are shown in the right panel. 
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Figure 16. The effect of photon sphtting on the final photon and pair spectra. The solid lines show the case where only photons with 
perpendicular polarizations are allowed to split (-L— >|| ||), and the dashed lines show the case where photons of both polarizations are 
allowed to split. The other cascade parameters are: Bp, 12 = 1000, Pq = 1 {a dipole field geometry is adopted), and 70 = 4 X 10^. The 
photon spectra are shown in the left panel and the pair spectra are shown in the right panel. For typical cascade parameters, the spectra 
when photon splitting is "turned off" in the simulation are nearly indistinguishable from the spectra when only -L— >|| || is allowed; we 
therefore do not present the turned-off case here (but see Fig. I12I I . 



rules for photon polarization modes, one-photon pair pro- 
duction into low Landau levels for the e^, and resonant in- 
verse Compton scattering from polar cap hot spots [with 
T = (0.3-3) X 10® K]. Both cascades initiated by a single 
electron (representing the entire cascade process described in 
Section [1]) and by a single photon (representing one branch 
of the cascade) are simulated, for a variety of initial ener- 
gies of the primary particle {■yom,^(? and eo)- We have made 
an effort to present our numerical results systematically, in- 
cluding empirical relations for the pair multiplicity (i.e., the 
number of electrons and positrons produced per primary 
particle) . 

Locally the cascade behaves very differently above and 
below the critical QED field strength Bq = 4.414 x 10^^ G 
(e.g., pair production is ~ 10 times more efficient for B > 
Bq but synchrotron emission is highly suppressed; see Sec- 
tion [4]). Globally, however, the cascade followed from the 
surface to the light cylinder behaves similarly regardless of 
surface field strength Bp. For example, we find that the to- 
tal number of pairs produced in electron-initiated cascades, 
as well as their energy spectrum, depends on the polar cap 
voltage BpP"^ but is weakly dependent on Bp alone. Addi- 
tionally, the total photon spectra for B > Bq and B < Bq 
have similar shapes over the energy range we consider in our 
simulation (10 keV-1 TeV); this is because curvature radia- 
tion, which dominates the photon spectrum over most of this 
energy range, is independent of magnetic field strength. Pho- 
ton splitting, a process that is only active for B > Bq , could 
potentially distinguish between neutron stars with high and 
low surface fields, by lowering the multiplicity of electrons 
and positrons Ne produced in the cascade to such a de- 
gree that the radio emi ssion mechanism can not function 
l|Baring fc Harding|200il l. However, we find that even if both 
photons polarized parallel to and perpendicular to the mag- 
netic field are allowed to split, photon splitting lowers A''^; by 
at most 50% for Bp,i2 = 1000 and by < 10% for Bp,i2 < 100. 
With the correct selection rule (-L— ^-H |1), the effect of photon 
splitting is even smaller. 



Our results show that a strongly non-dipole magnetic 
field with radius of curvature ~ 10® cm near the stellar sur- 
face (r less than a few stellar radii) can account for pulsars 
which lie below the standard death line for dipole fields (see 
Section [2| . Whether or not such a strongly-curved geome- 
try could form is another question altogether. It is doubtful 
that cascades initiated by resonant inverse Compton scat- 
tering in the open field region of the magnetosphere can ac- 
count for the "missing" pulsars, as the multiplicity of pairs 
produced in such cascades ( ^ 1 per primary electron) are 
too low for the current models of the pulsar radio emission. 
High-multiplicity cascades due to resonant ICS may still oc- 
cur along twisted lines in the closed field region (where the 
primary electrons never reach Lorentz fa ctors larger than 
70 ~ 10'^), as they do in mag nctars (e.g., 'Thompson et ah] 
I2OO2I: iBeloborodoy fc Thomps on 2007; Thompson 2Q08a ^ 
lBeloborodo^l2009l ). 



At any altitude in the magnetosphere, the photon spec- 
trum at low energies ( ^ 1 MeV) has a power-law shape 
with spectral index F — 2/3 (where dN/de oc (~^)\ this 
is true regardless of whether synchrotron radiation or cur- 
vature radiation dominates the spectrum, as both pro- 
cesses have the same low-energy shape. The photon spec- 
trum dN/de oc e"^/^ is harder than the hard X-ray spec- 
tra observed in several pulsars, which typically have pho- 
ton indices F ~ 1-2 (see, e.g., iKuiper fc Hermsen I I2OO9I : 
iBogdanov fc Grindlavl [2003 1 . One way to reconcile our re- 
sults with observations is to include an additional radiative 
process, cyclotron resonant absorption, in the si mulation. 
Such an approach is taken bv lHarding et al.l (|2005l ) (see also 
iHarding et al.l [20081 ) . who find that in addition to the syn- 
chrotron radiation emitted immediately upon creation, the 
particles emit an additional component of synchrotron 
radiation at large altitudes that dominates the low-energy 
spectrum. The efficiency of this high-altitude synchrotron 
emission is due to the large pitch angles of the , which they 
obtain through resonant absorption of radio photons that 
are beamed from some intermediate height in the magneto- 
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sphere. The high-altitude synchrotron emission will have a 
hard X-ray part with photon index Fio-iooo kcv generally 
different from 2/3, because the hard X-ray band lies near 
the peak of the emission (rather than in the low-energy tail 
as is t he case for the low-altitude emission). iHarding et al.l 
(|2005l ) find that for typical millisecond parameters the pho- 
ton index is in the observed range, 1 < Fio-iooo kov < 2. 

To fully incorporate the effects of cyclotron absorption 
into our simulation, we would need a model of the radio 
beam structure and spectrum, since the evolution in the 
pitch angle of each particle depends on the angle at which 
radio photons of the resonant frequency are incident on the 
particle an d their density. Empirical models of radio beams 
exist (e.g., Ran kinll 19931: lArzoumanian. Chernoff. fc Cordej 
l2002l : lKiiak^ Gil 200^) ; however, inclusion of such a model 
is beyond the scope of this paper. This means that we can 
not say anything in detail about the hard X-ray spectra 
of strongly-magnetized neutron stars, such as whether the 
high-altitude synchrotron comp onent dominates an d what 
its low-energy cutoff is (but see lHarding et al]|2005l for ex- 
amples of these spectra at several different cascade param- 
eters). Nevertheless, we can use the spectra from our 
simulation to estimate how the Fio-iooo kcV photon index 
varies as a function of the cascade parameters. For syn- 
chrotron emission from a distribution of particles, the 
hardness of the spectrum and the hardness of resulting 
photon spectrum are correlated [for a dN/ dE oc power- 
law distribution of the photon spectrum has an index 
given by the familiar expression F = (p -I- l)/2]. We find 
that the spectrum is harder, and therefore Fio-iooo kev 
is lower, for larger Bp or 70 or shorter P (i.e., anything 
that increases the energy of the cascade); additionally, for 
two neutron stars with the same polar cap voltage BpP~^, 
Fio-1000 keV is lower in the neutron star with the stronger 
surface field. Therefore, in pulsars where the cyclotron ab- 
sorption - synchrotron emission mechanism dominates the 
hard X-ray spectrum, the most-strongly magnetized pulsars 
will have the hardest spectra (lowest photon indices), all 
other parameters being equal. 

Luminous hard X-ray (from 20 keV to several hundreds 
of keV) emission has also be detected from a number of 
Anomalous X-ray Pulsars by INTEGRAL and RXTE (e.g., 
iKuiper et a"i]|2006l ). Possible mechanisms for this emission 
were discussed by Thompson & Beloborodov (2005) and Be- 
loborodov & Thompson (2007). Since the observed hard X- 
ray luminosity exceeds the spin-down power by several or- 
ders of magnitude, it must be fed by an alternative source 
of energy such as the dissipation of a superstrong magnetic 
field. We note, however, that the observed photon indices, 
F ~ 0.8— 1.8, are similar to (but slightly harder than) those 
of radio pulsars. Synchrotron radiation by secondary e^e~ 
pairs produced in a cascade similar to those studied in this 
paper could plausibly explain the observations. 

In constructing our simulations, we have attempted to 
rely as little as possible on the precise nature of the ac- 
celeration region. However, there are several key assump- 
tions that we have made that are only valid for inner 
gap accelerators (e.g., that the cascade begins at the neu- 
tron star surface). It has become apparent that models 
where particle acceleration occurs only in the inner gap 
can not account the observed high-energy gamma ray emis- 
sion from young pulsars. For example, the gamma-ray pulse 



profiles (e.g., the widely separated double peaks) of the 
six pulsars detected by EGRET already suggested that 
models of high-altitude gamma-ray emission were required. 
More recently, observations of the Crab pulsar by the Ma- 
jor Atmospheric Gamma-ray Imaging Cherenkov Telescope 
(MAGIC) and around 50 pulsars with "above average" spin- 
down powers by the Fermi Gamma-ray Space Telescope 
have revealed that the high-energy tails of the photon spec- 
tra fall off exponentially or more s lowl y than exponential 
(iAliu et al.ll2008l : lAbdo et al.ll2009allbll3 . l2010h . while inner 
gap models predict that the tails fall off faster than ex- 
ponential (see, e.g., Fig. I12p . The outer gap model and 
the high-altitude version of the slot gap model have been 
successful in explain i ng the 7-ray pulsar light curve s (e.g. , 
IWatters et ahl l2009l: IVenter. Harding, fc GuillemotI l2009l : 



but see lBai fc Spitkovskvll2009l ')? 

Although an outer gap or slot gap accelerator model 
is required to explain gamma-ray observations, our results 
based on the inner gap model still have general applica- 
bility. First, even when the acceleration region is located 
in the outer magnetosphere a significant fraction of the 
pair creation mus t occur in the inner magnetosphere (e.g., 
ICheng et al.ll2000l ). Indeed, pulsar radio emission is strongly 
constrained to originate from well inside the light cylinder 
radius; the only way to generate highly coherent radio emis- 
sion is to have copious plasma produced by vigorous pair 
cascades (e.g. , .Melrose 2004 ; .Lvubarskv..2008.1 . Second, both 
inner gap and outer gap accelerators may exist in a neutron 
star simulta neously, wheth er as one extended acceleration 
region (e.g., iHirotanil I2OO6I ) or on different field lines. Ra- 
dio observations suggest that radio emission can come from 
intermediate field lines, neither along t he pole of the sta r 
nor at the edge of the open field region (|Gangadharall2004l '): 
if this is true, it provides further support for an inner gap 
origin of the plasma, since the plasma generated by an 
outer gap model is formed on field lines close to the last open 
field line. In addition, as discussed above, hard X-ray (10- 
1000 keV) observations of both pulsars and magnetars sup- 
port a magnetosphere model where the hard X-ray emission 
is dominated by radiation from the e^e^ pairs generated 
by inner gap cascades. Third, it is unlikely that the death 
lines for the inner gap and outer gap mechanisms overlap 
completely. There should therefore be regions of the P-P di- 
agram where the outer gap mechanism is excluded but the 
inner gap mechanism still functions; these regions, if they 
exist, will be located near the inner gap death line, where 
pulsars have very low spindown powers. Although observa- 
tions show that the gamma-ray spectra is dominated by slot 
or outer gap emission for pulsars with moderate to large 
spindown power, they have not yet ruled out an inner gap 
origin for pulsars with low spindown power. 
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APPENDIX A: RESONANT INVERSE 
COMPTON SCATTERING 

Here we calculate the photon scattering rate for the reso- 
nant inverse Compton process, using the simplified model 
of an electron positioned directly above the magnetic pole 
(hot spot) and traveling radially outward; see Fig. lAll The 
resonant cross section for inverse Compton scattering, in the 
rest frame of the electron before scattering, is 



(Al) 



where e' = 7ei(l — ficostp), is the photon energy in the 
"lab" frame, and is the incident angle of the photon with 
respect to the electron's trajectory. This cross section is ap- 
propriate even for B > Bq, since the resonant condition 
e' = Ec holds regardless of field strength [though at the 
highest field strengths the prefactor e^h / nif.c in Eq. (lAlf) 
serves only as an up per bound to the exact expression; see 
iGonthier et al.ir2000l ]. The polar hot spot has a angular size 
of ^spot; this sets a maximum value for ^ of 



cos t/jc; 



iicOS Ospot 



^r^ + -2rRcosespo 



(A2) 



where r is altitude of the electron (from the center of the 
star). The intensity of emission from the hot spot is 



.(r) = ^ 



Ip <tpc 



otherwise. 



Therefore, the photon scattering rate for the resonant ICS 
process above a polar hot spot is given by [see, e.g., Eq. (B3) 
of ML07] 



dN, 



ph 



dt 



J dQi J dti (1 — /3C0S?/)) ^ — ^ CJre 



(A4) 



dQi 



kT 



B, 



!»=fc/[7(l-/3cosV)] 



Pq In 



(A5) 

1 _ e-ec/[7(l-/3cos^„it)'=T] ' 

(A6) 



APPENDIX B: PAIR PRODUCTION 
Bl Kinematics 

Consider the pair production of a photon with energy e and 
angle ij) (the photon - magnetic field intersection angle). In 
the frame where the photon is traveling perpendicular to the 
local magnetic field direction (this "perpendicular" frame 
moves at the velocity ccos-i/) relative to the "lab" frame), 
the photon has energy e = esimp, and energy conservation 
demands 



(Bl) 



where Ej and E'^. are the energies of the electron and the 
positron and is the momentum along the magnetic field 



electron 



photons 



pole 




J 



Figure Al. A simplified picture of the ICS effect on the electron. 
The electron is assumed to be directly above the magnetic pole 
and travel radially outward. 



of either particle (p^ 



From this we find 



1 - (j + k)fiQ + (i - k)' 



where x = t /{2m,^c?), and 



E' 



TUeC 



= m^c\lx^ + {k- j)PQ + U - fc)2^ 



(B2) 



(B3) 



(B4) 



In the "lab" frame, the energies of the electron and positron 
at the moment of pair creation are given by 

E = ^^(E' ±p',c cos -iP). (B5) 
sm ?/) 

In our simulation one particle is randomly assigned the '4-' 
energy and the other the ' — ' energy, with equal probability 
of either outcome. 



B2 Photon attenuation coefficients and optical 
depth 

In the perpendicular frame, the first three attenuation coeffi- 
cients for II and the first two non-zero attenuatio n coefficients 
for _L photons are (|Daughertv fc Hardindl 19831 ) 

1 Pq 



R\\,oo — 



R'u 



11,01 



2 X 



2ao xVa;^ - 1 
2 + /3q 



x> 1; 



(B6) 



2ao 



Rlnn=2X 



1 2x 



^2 - 1 - /3q + 

X > (l + v/1 + 2/3q) /2 ; (B7) 
1 + , 



2x^ 



2«o /3q J,2_^_2Pq + 'A 



X > 1-f \/1 + 4/3q /2; (B8) 
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R'±,oi = 2 X 



2ao 2x^ 



2x' 



+ 



x> (1 + V1 + 2/3q)/2; (B9) 



R\ 



2 X 



2ao 



a; > (l + \/1 + 4/3q) /2 . (BIO) 

Note that -Rx.oo — 0- I n the above expressions, th e at- 
tenuation coefficients of iDaughertv fc Harding! (|l983l ) are 
multiplied by a factor of two for all channels but (00) 
[i.e., in Eqs. HB7|) - HB1Q|| ]. since we are using the convention 



r,/ r)/(DH83) , r>/ (DH83) 



2R 



J (DH83) 



,fc T ..fe^ - -..^.fc for j / k. 
We now examine the conditions for pair production by a 
1 1 -polarized photon; the analysis is similar for a _L-polarized 
photon and yields the same result. The optical depth for 
pair production is 



dsR{s) = 



"'0 



ds R'{s) sin jp . 



(Bll) 



We assume tl) <^ 1, which is valid since most photons that 
can pair produce will do so long before the angle t/j ap- 
proaches unity (only photons with energies e ~ 2meC^ in 
the lab frame must wait until i/; ~ 1 to pair produce). In 
this limit, we have sin?/; ~ s/TZc, so that x and s are related 
by 

X ~ ^— . (B12) 

Tic 2mec2 ^ ^ 

Let soo to be distance traveled by the photon to reach the 

first threshold x — xqo = 1, and soi to be the distance 

traveled by the photon to reach the second threshold x = 

xoi = ^1 -I- ^/l + 2/3q^ /2. The optical depth to reach the 

second threshold for pair production is 



roi = 



dsii||,oo(s) 



2ao 



2meC 



dx 



•c\fx 



(B13) 



(B14) 



= 9.87 X 10" 



11 / 7^c \ ( 100 MeV 



V 10» cm ) \ 



T{Pq), (B15) 



where 



T{Pq) = /3q 



dx 



x^fxP-'- 



(B16) 



We plot Toi as a function of magnetic field strength in 
Fig. [ml for e = 100 MeV and R^ = 10* cm. 

From Fig. [m] we see that pair production occurs in the 
(jfc) = (00) channel (tqi > 1) when 

B > B„it ==; 3 X 10^^ G. (B17) 

Because of the steep dependence of t on _B for _B ~ 3 x 
10^^ G, the value of Bcrit does not change much for different 
parameters e and Tt^- For example, -Bcrit = 3 x 10^^ G for 
e = 100 MeV and R^ = 10* cm, and Bcrit = 7 x 10^^ G 
for e = 10* MeV and i?c = 10^ cm. Figure [Bl] also shows 
that for _B< Bcrit, the optical depth roi is much less than 
unity. The same result is found for the optical depth from the 
second to the third threshold, and for higher thresholds. The 




0.2 0.4 0.6 0.8 1 

Pq=B/Bq 

Figure Bl. A plot of the function T(/3q) as given by Eq. I IB16II 

(right axis) and the optical depth to reach the second threshold 
for pair production, tqi (left axis), as a function of for e = 
100 MeV and Tic = 10* cm. 



pair production process can therefore be divided into two 
regimes: for B<i?crit, photons must travel large distances 
before pair producing, at which point the resulting pairs will 
be in high Landau levels; for B > Bcrit, photons pair produce 
almost immediately upon reaching the first threshold, so 
that the pairs will be in low Landau levels (n<2). 



APPENDIX C: EMPIRICAL RELATIONS FOR 
THE NUMERICAL RESULTS 

In this section we justify several of the empirical relations 
given in Sections 14.11 and 14.21 In the derivations below we 
treat the radius of the emission point for the primary pho- 
ton, re = r'o,ph/(10^ cm), as a free parameter. This allows 
the results to be applicable both to Section 14.11 where we 
assume that re, = 1, and to Section [4.21 where re > 1. 



CI Photon-initiated cascades 

At low fields (B< Bcrit 3 x 10^^ G), we can use Eq. ^ 
for the attenuation coefficient, which implies that pair pro- 
duction occurs when xfiq ~ 1/15-1/10. In the following 
we will use x^q — 1/10, appropriate for photon energy 
e ~ 10'' MeV. To initiate an effective cascade, a photon must 
pair produce before traveling a distance Sph,max ~ 0.5ro,ph 
[Eq. (|31|l ]. Therefore the minimum photon energy for cas- 
cade is given by 

emin O.Srp^ph ^ 1 ^^^^ 



2meC^ Tic 



WI3q 



where /3q and TZc should be evaluated at the pair creation 
point r ~ Sph.max 4- ro,pii = 1.5ro,pi,. Since /3q = B/Bq 
0.02Bp,i2(r/ii)~^ and 7^c ~ Kflls^JrjR cm, where Bp,i2 
and TZs, are the magnetic field strength (in units of lO^'^ G) 
and curvature radius (in units of 10* cm) at the surface, we 
find 



emin ~ 3000B-^27^8r6^^ MeV 



(C2) 



For a photon injected from the surface (re = 1), this reduces 
to Eq. (gSJ. 

In each generation of pair production, a photon of en- 
ergy e and angle ?/) creates an pair, each with energy 
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7meC = 0.5e and traveling in the same direction as the 
photon. The electron and positron radiate synchrotron pho- 
tons until they reach an energy of 7||meC'^ = met? j sin-!/) = 
^meC? /x ~ 0.1-Bi2e, since x = "ysmip ~ l/(10/?<3). The 
synchrotron photons have a characteristic energy esR = 



1.57 {s\n%li)PQmeC 



0.075e. Therefore, in each genera- 



tion, the energy of the photons and pairs drops by a factor 
~ 0. 075 w hile the number of pa rticles increases by 1/0.075 
fsee Imbschman fc AronsI l2001ai l. When RIGS is inactive 
(e.g., when T < 10* K), such that the only mechanism for 
energy loss is synchrotron radiation, the final pair spectrum 
consists of an pair, each with energy 

~ 0.1i3p,i2rg ^eo , (C3) 

created by the primary photon, and a power-law distribution 
of pairs with index p ~ 2 (where dN/dE oc E~^) extend- 
ing from ~ 0.075i?max down to ~ 0.1-Bi2emin (see Fig. (B]). 
Hence the total energy of the cascade pairs produced by the 
primary photon (of energy eo) is 

'0.075eo' 



Stot ~ 2_Emax + 0.1i3p,12rg eminTlE lu 



(C4) 



which reduces to Eq. (|52|l for re = 1. 

For the range of high fields considered in this work 
(Scrit < 1000), pair production occurs when x = 1 (for 
X = xoo) or X ~ 1-4 (for x = xoi; see Section [3.2|l . For a 
photon emitted from the point (ro.ph, ^o,ph), the maximum 
intersection angle between the photon and the magnetic field 
is given by [Eqs. HHI, (HZI, and (^1]) ] 

siUT/jmax ^XlX (Slhph)] - X {So,ph) 



(C5) 



arctan 



■ tan 



X + arctan 



tan 9o. 



ph 



(C6) 



for "small" angles &o,ph < 0.6, 
sin^ 

max 

< sin 6'o,ph • (C7) 

For 9 < 0.6 the local radius of curvature is given by TZc — 
1.3r/ sin 61; therefore, using x — esinip/{2me?) ~ 2 we have 

1/2 



MeV. 



(C8) 



e,nin ~ 2007^8r■g 

which at re = 1 is Eq. (f54|) . 

At high fields, each electron-positron pair is created 
with energies ■ynrieC^ ~ 0.5e (for || photons) or ~ Q.bt/xQi 
and ~ 0.5e(2 — l/xo\) (for _L photons). After synchrotron 
radiation the electron and positron have energies 0.5e or 
O.Se/soi. Therefore, the final energies of the electron and 
positron created by the primary photon are given by 
Co eo 



E„ 



or 

2x01 2x00 



• O.leo or 0.5eo • 



(C9) 



In order for resonant ICS to modify the cascade spec- 
tra, the electron (positron) must have an energy (after syn- 
chrotron radiation) larger than the minimum energy at 
which RICS is effective, i.e., -Emax > 0.37critmeC^ [where 
7crit = ec/kT; see Eq. (|53|) ]. which implies 



eo > JOBpaaTg-i MeV. 



C2 Electron-initiated cascades 



(CIO) 



The spectrum of curvature photons extends from approxi- 
mately one photon at 

emax ~ lOecR (70) = 3 X 10S?7^8 ' MeV , (CU) 



where 77 is the primary electron's initial Lorentz factor 70 

in units of 10 up to a maximum of ~ 6 x lO^T^-g^^ photons 
1/2 

at epcak ~ 67?.g MeV. For photon energies below tpoak the 
spectrum is a power law with spectral index F = 2/3 (where 
dN/dE oc E~^), characteristic of curvature and synchrotron 
radiation at low energies re.g.. lErbeilll966l ). 

At low fields B < -Bcrit the pair spectrum extends from 

E 

max ,12£max 3007?Bp,l27^8 ' MeV (C12) 



down to 



E^in ~ O.lBizemin ~ 3007^8r6-'/^ MeV ; (C13) 
at high fields 

£;n,ax 0.5e,nax ~ 6 X 10^7?7^^' MeV (C14) 

and 

£;n.in 0.5en.in/a:oi ~ 2Qnsr-^^^ MeV , (C15) 
(see Section [Cl|l . 
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